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TEST SECTIONS FOR SMALL THEORETICAL 
WIND-TUNNEL-BOUNDARY INTERFERENCE 
ON V/STOL MODELS 

By Ray H. Wright 
Langley Research Center 

SIMMARY 

A wind-tunnel test section with closed upper wall, slotted side walls, and open 
lower boundary was found theoretically to produce zero tunnel-boundary lift interference 
on a small wing with horizontal wake mounted at the center of the test section. For this 
test section the variation of the interference with angle of the vortex wake behind a high- 
lift- coefficient model was not large. Because of the small slot widths required for zero 
interference and of the effects of boundary layer, the theory is regarded as unreliable for 
predicting the slot widths required; however, the variation of the interference with the 
slot width for widths somewhat greater than those needed for zero interference was found 
to be small. The interference in the region likely to be occupied by the tail of a model 
was investigated in some detail and was found to change with slot width and with wake 
angle more strongly than did the interference at the lifting element. A limited investiga- 
tion of the lift interference in a test section with closed upper wall and slotted side and 
lower boundaries was made to obtain the theoretically indicated slot width required for 
zero lift interference at a center- mounted wing with the wake horizontal. 

INTRODUCTION 

Wind-tunnel test-section configurations designed to minimize the boundary inter- 
ference on the flow over test models have been in use for many years. (See, for example, 
ref. 1.) For high-speed models the interference of greatest concern is that due to solid 
blockage, because this type of interference increases strongly as the test Mach number 
increases toward unity. On the other hand, for testing models of high- lift- coefficient 
vehicles, such as vertical-take-off and vertical- landing craft (VTOL) and short-take-off- 
distance and short- landing- distance craft (STOL), which generally operate at relatively 
low Mach numbers, the lift interference is of first importance. With these high-lift- 
coefficient models the problem is complicated by the fact that the wake associated with 
the lift may vary from approximately line vortices extending almost horizontally down- 
stream from the model to a jet extending vertically below, or even upstream from, the 


model. The wind- tunnel-boundary interference associated with idealizations of heli- 
copter wakes in the presence of several configurations of open and of closed vertical and 
horizontal boundaries has been treated in references 2 and 3. 

Some insight into the lift interference with high- lift- coefficient test models can be 
obtained by examining the calculated interference correction factors for helicopter rotors 
presented in reference 3. For most test- section configurations, the interference varies 
widely, depending on the angle made by the wake of vorticity with the vertical. However, 
as may be seen from figures 18 to 21 of reference 3, if only an open lower boundary is 
present, this variation (with the model some distance from the open boundary) is small. 
The calculated interference due to lift is greatest for the wake extending vertically below 
the rotor; however, it is believed that for this vertical wake, the calculations of refer- 
ence 3 are unreliable, because the wake, which for this condition is in the form of a jet, 
might be expected to break through the open boundary of the stream rather than to follow 
along its edge as assumed. In a wind tunnel with boundaries not completely open, as the 
jet approaches the vertical and the tunnel stream velocity approaches zero, the interfer- 
ence velocity at the rotor must be an upwash (for positive lift) because any solid boundary 
restricts the flow. The open- lower-boundary condition can be approximated with a many- 
slotted bottom wall, and the multiple- slotted bottom wall might be preferable to the com- 
pletely open lower boundary in order to stabilize the flow. 

The purpose of the present investigation is to extend the small-disturbance theory 
of wind-tunnel lift interference on high- lift- coefficient models, such as reported in refer- 
ences 2 and 3, to some test- section configurations that might be expected to minimize the 
theoretical lift interference for the general class of high- lift- coefficient models over 
their operating ranges. The configurations investigated are of the mixed (open and 
closed) boundary type. First to be treated is a test section with top and side walls 
closed and lower boundary open. The model, assumed to be mounted on the vertical cen- 
ter line of the test section, is represented by a vortex doublet trailing horizontally down- 
stream. A vortex doublet is the limit as two rectilinear parallel line vortices of equal- 
magnitude but opposite- sense circulation approach each other, the product of circulation 
with their distance apart remaining constant. This product is called the strength of the 
doublet. Representation of the wing by a trailing- vortex doublet omits the bound part of 
the vortex system, but the interference at the lifting element is not affected by this 
omission. The results are therefore valid for a small, lightly loaded wing. Calculations 
of the interference are made for rectangular test sections having height-width ratios of 1 
and 1.5. 

The second configuration treated is like the first except that the side walls are 
slotted with four equispaced slots. Added flexibility is thereby available for obtaining 
desirable interference characteristics. 
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In order to remove the limitation of small lift loading, this second test- section con- 
figuration with height-width ratio of 1.5 and slot width corresponding to small downwash 
interference at the model is investigated with a helicopter rotor mounted at the center. 

As in references 2 and 3, the rotor wake is represented by a source- sink doublet dis- 
tribution. This representation again implies the small-model assumption. Even with a 
^ small model, the small -disturbance assumption is violated over part of the open lower 
boundary, because the now- no- longer- horizontal wake approaches and even penetrates 
that boundary. The validity of the theory may be expected to be degraded progres- 
• sively with expansion of the region of the boundary over which the small-disturbance 
assumption is seriously violated. The small-model assumption therefore takes on added 
importance. On the other hand, provided the small-model assumption is satisfied, the 
lift-interference theory for the helicopter applies also for other high- lift- coefficient 
models, and for a horizontal wake it is exactly equivalent to the interference theory for 
a small, lightly loaded wing. Moreover, if the wake is deflected slightly below the hori- 
zontal, it intersects the open lower boundary so far from the model that the violation of 
the assumed boundary conditions near the intersection region can have little effect on the 
interference in the vicinity of the model. The helicopter interference theory should there- 
fore be valid for calculating the lift interference on general small high- lift- coefficient 
models with wake deflected not too much from the horizontal. The wake deflection can 
be measured as part of the tests or it may be estimated from the lift, the area and geom- 
etry of the rotor or lifting surface, and the tunnel stream velocity. (For a short discus- 
sion of this problem, see ref. 4.) Unfortunately, although a separate solution is obtained 
for an exactly horizontal wake, the general solution of the helicopter rotor interference 
with deflected wake involves nonuniformly convergent integrals; thus, in the calculations 
the convergence problem becomes increasingly severe as the wake approaches the 
horizontal. 

As an effort toward completeness and toward provision of some theoretical results 
for comparison with experiment, several other test- section configurations are considered 
for which the theoretical lift interference on wings (small deflection of the wake) may be 
reduced to zero. These configurations include the interference of the model for small 
wings centered over solid lower boundaries with top and side boundaries open and height- 
width ratios of 0.5, 0.66, 1.0, and 1.5. The possibility of reducing to zero the lift inter- 
ference of small wings by means of slots in all four walls is also mentioned and appro- 
priate references are given. Finally, the lift interference is calculated for a lightly 
loaded center- mounted wing spanning 0.58 of the width of a test section of height-width 
ratio 1.432, with closed top and with sides and bottom slotted with four equispaced slots 
each. The wing is represented by line vortices trailing downstream from the tips. The 
calculations are made for a range of slot widths including that width for which the inter- 
ference of the model is theoretically zero. 
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Of these several theoretical developments only that for the helicopter rotor is 
applicable for calculating the interference at the tail. Several studies are therefore made 
in order to understand the various contributions to the upwash interference at the tail and 
to obtain a qualitative estimation of that interference for the test section with solid top 
and slotted sides and bottom. These studies include consideration of the tail interference 
corresponding to the trailing vortices, investigation of the effect of below- center position 
of the trailing vortices in the closed-top slotted- side-wall and slotted-bottom configura- • 
tion, and qualitative consideration of the tail interference of solid top and open bottom and 
of slotted side walls in the presence of the bound vortex. Finally, the helicopter theory 
is adapted to calculate first, the effects of infinite parallel slotted side walls and second, 
the effects of the closed top, open bottom, and slotted side-wall configuration on the 
upwash interference downstream from the model. For the first of these the calculations 
are made for walls slotted to produce zero interference of the lifting element. For the 
second the calculations are made for two configurations of slotted walls, one for zero 
interference at the lifting element and the other corresponding to a small downwash 
interference. 

In addition to .the restrictions imposed by small-model and small-disturbance 
assumptions, unavoidable uncertainties in boundary conditions limit the applicability of 
the theoretical treatment such as the possibility, with very high lift, of separation from 
the upper wall. Because of such uncertainties, experiments are required to check the 
theory and to determine the desirable slot widths. The theory should nevertheless fur- 
nish useful guidance for experimental investigations directed toward developing test sec- 
tions to reduce the tunnel-boundary interference on high- lift-coefficient models. Since 
accurate corrections for the interference do not seem likely, such mixed-boundary test 
sections are believed to hold little promise for increasing the model size that can be 
satisfactorily tested in a test section of given cross section. (For some discussion of 
the model- size problem, see refs. 4 to 7.) The more likely benefits are the elimination 
of all corrections with properly small models and possibly an improvement in the uni- 
formity of the effective test flow in the region of the models. 

SYMBOLS 

Symbols pertinent to text and figures are given in this list. Appendixes A to D 
have separate symbol lists. 

A area of helicopter rotor disk 

A vector area of rotor disk 

b width of test section 
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c 

c L 

g 

h 

l 

- P 

r 0 

r 

S 

s 

u 

V 

V 

w 0 


cross-sectional area of test section 

lift coefficient of test model, positive for upward lift 

distance from model position to lower boundary of test section 

height of test section 

restriction constant (defined by eq. (B2)) 

field point 

ratio of slot width to distance between centers of panels of a slotted wall 
vector distance from element of surface enclosed by vortex ring 
model area on which Cl is based 

distance from rotor disk to vortex ring 

streamwise interference velocity, positive in tunnel stream direction 
tunnel stream velocity 

upwash interference velocity, positive upward 

velocity at rotor disk induced by the lifting rotor (taken positive for positive 
lift) 


X,Y,Z 

x ,y>z 

a 

r 


5 R,x 

6 R,z 


0 


rectangular Cartesian coordinate axes 

rectangular Cartesian coordinates 

angle of attack of helicopter rotor-tip-path plane 

circulation, positive for positive lift 

upwash interference factor for wings 

streamwise interference factor for helicopter rotor 

upwash interference factor for helicopter rotor 


5 


strength of vortex doublet 


4 

p density of test medium 

X angle between axis of helicopter rotor disk and axis of cylindrical vortex 

sheet extending from rotor, designated skew angle (positive 
counterclockwise) 


TEST SECTIONS 

Top and Side Walls Closed; Lower Boundary Open; 

Horizontal Trailing Vortex Doublet 

The boundary conditions for closed top and side walls and open lower boundary are 
simply satisfied by means of images as in reference 8. Figure 1 indicates the required 
arrangement of line doublets extending downstream from the position of the test model. 
The original lifting doublet at the position of the model is included in the figure. In fig- 
ure 1, b is the width of the test section, h is its height, and g is the distance of the 
model from the open lower boundary. The model is located on the vertical center line 
of the tunnel. An equation for calculation of the boundary- induced upwash velocity v at 
the position of the model is developed in appendix A. In this and other equations devel- 
oped in this investigation, no account has been taken of compressibility effects because 
in tests of V/STOL models, the speed of the test stream is normally so small that such 
effects would be negligible. 

It is convenient to express the interference in terms of the interference upwash 
factor where 6 W is related to the upwash interference velocity v by 

tanAas ^ B! "*r 8 w (1) 

where 

v upwash velocity, positive upward 

V tunnel stream velocity 

C cross-sectional area of test section 

Cl lift coefficient of model, positive for upward lift 

S area on which Cl is based 

A a effective change of angle of attack due to upwash velocity 
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2-(2h) - - 


+ + 4 + 4 + 4 - - - 

2-(2h)-2g - - - 4 4 + + 4 4 4 - - - 

I 

l 

2h - • - 
2h-2g 


0 + 4 

-3b -2b 

-2g 4 4 



-2h - - - - - - 

l 

-2h-2g - - - 


-2-(2h) - - - 4 4 4 

-2(2h)-2g 4 4 4 



4 4- 


4 4 


Figure 1.- Arrangement of vortex doublets satisfying boundary conditions for tunnel closed at top and on sides 
and open on bottom. +, doublet with upward lift; doublet with downward lift. 
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In order to express 6^ in terms of the doublet strength 4 instead of in terms of Cl, 
note that 

VC L S = 2 4 

and 


C = hb 

so that equation (1) gives 

«. hb 

5 » = 2^ v (2) 
Use of equation (A14) in equation (2) shows that the interference upwash factor 6^ is a 
function of the test-section height-width ratio h/b and of the ratio of model height above 
the lower boundary to the tunnel width g/b. 


The interference upwash factor 6^ has been computed by use of equations (2) 
and (A14) for the square test section ^ = 1 and for a rectangular test section with 
j-j = 1.5 with several values of g/b. The computed values of 6w are shown plotted 
against g/b in figure 2. 


Figure 2 shows that with the model near the center of the tunnel, the interference 
&W is approximately the same as that in the closed tunnel, that is, 6w = 0.137 for the 
square test section. On the other hand, if both top and bottom boundaries had been open, 



Figure 2.- Interference upwash factor for vortex doublet in wind tunnel with top and side walls closed 
and bottom open. Vortex doublet on vertical center line; wake horizontal. 
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the interference would have been 6w ~ -0.137 for the square test section, as can be 
seen from figure 8(a) of reference 9. This effect of opening the top as well as the bottom 
is due to the fact that with this configuration, the whole tunnel stream is deflected down- 
ward. This behavior shows the importance of preventing separation from the upper wall, 
since the interference factor might then approach the open-tunnel value. 

Figure 2 shows that in order to obtain zero interference in a test section with 
closed top and side walls and open lower boundary, the model must be placed relatively 
close to the open boundary in a region where the interference varies strongly with model 
position. In addition, it is to be expected that the interference varies strongly in the 
region surrounding a model placed in this position. This situation is unsatisfactory. 

Since the interference at the center is an upwash and since experience indicates that slots 
produce an interference effect of the opposite sign, slotted side walls to reduce the inter- 
ference to zero for the model mounted at the center of the test section are suggested. 


Top Wall Closed; Side Walls Slotted; Lower Boundary Open; 

Horizontal Trailing Vortex Doublet 

A lifting model of small span mounted on the vertical center line of the rectangular 
test section is again replaced for the purpose of calculating the boundary- induced inter- 
ference by a vortex doublet trailing horizontally downstream. The top is again solid and 
the bottom open, but the side walls are slotted with four equispaced slots so arranged that 
the center of a (solid) panel 
intersects the horizontal center 
line and a half-panel adjoins 
each corner. (See sketches in 
figs. 3 and 4.) The top and bot- 
tom boundary conditions are 
again satisfied by images 
arranged above and below the 
model as in figure 1, but the 
side-wall-slotted boundary con- 
dition cannot be satisfied by the 
b- interval repetition of vortex 
images of the same sign, which 
was effective for satisfying the 
solid-boundary condition. 

Instead* the variable in the ver- 
tical direction is eliminated by 

Figure 3.- Upwash interference factor 6w for vortex doublet in wind tunnel 
taking the exponential integral with top wall closed, bottom wall open, and side walls slotted, each with four 

equispaced slots. Vortex doublet on vertical center line*, h/b = l ; wake 
horizontal. 
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Fourier transform of the original disturbance potential, of the governing Laplace equa- 
tion, and of the average (homogeneous) boundary condition. Solution of the transformed 
system and inversion of the solution give the potential function for the doublet in the 
presence of the slotted walls. The total interference potential is obtained by summation, 
and from this potential, the upwash interference factor 6^. is derived. The complete 
development is given in appendix B. 

The upwash interference factor 6 ^ at the position of the lifting model has been 
computed from equation (B15) for the square test section and for a rectangular test sec- 
tion with £ = 1.5 for three heights g/h of the model above the lower boundary and is 
b 

shown in figures 3 and 4 as a function of the ratio of open to closed area of the slotted 
sides Tq. The slotted side walls are seen to be effective in reducing the upwash inter- 
ference to zero for a model located near the center of the test section; that is, for 
£ = Moreover, in the range of Tq values corresponding to &«, = 0, the interference 
for the centered model is seen to be relatively insensitive to slot width. These charac- 
teristics for the small model with vortex wake trailing horizontally downstream are 
regarded as favorable. Next to be investigated is the behavior in this test section of a 
high- lift model, a helicopter rotor, for which the vortex wake may not trail approximately 
horizontally downstream. 



Ratio of slot width to distance between centers of panels, r Q 

Figure 4.- Upwash interference factor 6w for vortex doublet in wind tunnel with top wall closed, bottom wall 
open, and side walls slotted, each with four equispaced slots. Vortex doublet on vertical center line*, 
h/b = 1.5; wake horizontal. 
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Top Wall Closed; Side Walls Slotted; Lower Boundary Open; 

Helicopter at Center of Test Section 

The helicopter rotor is assumed to be mounted at the center of the test section and 
as in references 2, 10, and 11, its disturbance is represented by a straight elliptic cylin- 
drical vortex sheet extending from the rotor. The rotor disk makes an angle of attack a 
with the oncoming wind-tunnel stream velocity V as indicated in figure 5, and the vortex 
cylinder is swept downstream by the skew angle y; thus, the axis of the vortex cylinder 
- makes an angle y - a with the downward- drawn vertical. Note that y is here mea- 
sured from the downward-drawn normal to the rotor disk rather than from the negative 
Z-axis as in references 2 and 3. As in the preceding treatment of the wing, the model 
is assumed to be small. This assumption permits the simplification of replacing the 
vortex cylinder by a source- sink doublet line. The calculations of reference 3 show that 
as the vortex wake is swept back toward the horizontal, that is, y - a = 90°, the boundary- 
induced upwash with the helicopter rotor approaches that for the wing producing the same 
lift. It is therefore unnecessary to investigate the interference on the helicopter rotor 
for values of y - a approaching 90°. This investigation is therefore mainly confined to 
values of the angle for which the vortex wake intersects the lower boundary within the 
test section. The vortex wake, which now resembles a jet, is assumed to break through 



Figure 5.- Elliptic cylindrical vortex sheet representing disturbance^ due to helicopter rotor. Vector area of vortex ring is 

taken equal to vector area A of rotor disk. 
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the lower (open) boundary and to have no further influence on the flow about the test 
model. This assumption seems reasonable if the dimensions of the open space sur- 
rounding the test section are large in comparison with the diameter of the rotor. It is 
not, of course, correct if the wake intersects the closed wall of the diffuser, but the 
equations herein derived should, even in this case, give a reasonable approximation to 
the boundary- induced interference provided the diffuser entrance is sufficiently far down- 
stream, as it should be, to exert small influence on the flow about the test model. 

The boundary conditions assumed are essentially those of the preceding section, 
namely zero velocity normal to a solid wall and zero pressure increment at an open 
boundary. Even for a wing the boundary condition at an open boundary is uncertain, 
largely because the position of the boundary is unknown. For a helicopter with wake 
approaching the open lower boundary, it is obviously far less certain, and with the tunnel 
stream velocity V approaching zero (and x — 0)> the assumed boundary condition at the 
open lower boundary becomes completely invalid. The lift- interference theory for the 
helicopter is therefore generally less reliable than that for the small, lightly loaded wing. 
Nevertheless, it should furnish useful indications, even for the x — 0 condition, pro- 
vided the helicopter rotor is not too close to the open lower boundary. 


Equations for computing the interference factors 6 R and 6 R z for the heli- 
copter rotor in the test section with closed top, slotted side walls, and open lower bound- 
ary are derived in appendix C. Here 


and 


k - u 
°R,x " w f 


C 

A 


(3) 


fi = 

°R,z Wq A 

where 


(4) 


u boundary -induced interference velocity in tunnel stream direction 

v boundary- induced interference velocity in vertical direction, positive upward 

Wq average velocity induced by lifting rotor at rotor disk 

C cross-sectional area of test section 

A area of rotor disk 

In this paper the interference velocity components u and v are derived from a single 
potential rather than being considered separately as in references 2 and 3. Thus the 
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drag terms of these references are already included. (Note that wq is the resultant 
induced velocity at the rotor disk and at a = 0 it is the negative of the corresponding 
quantity used in refs. 2 and 3.) If the helicopter is operating at considerable forward 
velocity, its lift is approximately 2pVwgA whereas the lift of the wing is ^pV^ClS. If 
for this operating condition the ratio of upwash interference velocity v to total lift is the 
same as that of the wing, it therefore follows from equations (1) and (4) that 

6 R,z = 45 w 

In the expression for the lift, p is the density of the test medium. 

Interference factors for the helicopter rotor at the center of the test section with 
closed top, slotted side walls, and open lower boundary and with height-width ratio h/b 
equal to 1.5 have been computed from the equations of appendix C by use of an IBM 7094 
electronic data processing system. For this calculation a homogeneous boundary condi- 
tion on the side walls corresponding to the four equispaced slots in each wall with an open 
ratio Tq of 0.0455 was assumed. Because of a mistake in calculations, this value of 
Tq was originally believed to correspond to zero interference on the small wing model 
mounted at the center of the test section (this incorrect result has been quoted in ref. 4); 
but as seen from figure 4, it actually corresponds to by? = -0.026. The calculation was 
made for angles of attack a of -10°, 0°, 10°, and 20° and for several values of the skew 
angle y. The streamwise interference factor 6 n was computed at the position of the 

test model £ = 0 only. The upwash interference factor 6 D was computed at the 
h k,z 

model position ^ = 0 and at the downstream positions ^ = 0.1, 0.3, and 0.5. The cal- 
culated interference factors are shown in figure 6 as functions of the angle y ~ a between 
the wake axis and the downward-drawn vertical. The values are, in general, not com- 
pletely negligible, but every curve approaches or crosses zero for some value of y - a 
and the variation does not seem to be unreasonable. For instance, with a = 0, ^ = 0, 
an upwash interference 6-r, -0.104 occurs at y = 0. Such an upwash interference is 

required, because any boundary other than a completely open boundary places a constraint 
on the flow for this hovering mode of operation. As y increases, the value of St, at 
- = 0 decreases, becomes negative, and then appears to be tending back toward zero as 
it should if for values of y - a approaching 90° the interference factor on the helicopter 
rotor is to approach 6,, =-0.104 corresponding to ^ = -0.026 for the wing. It is 

x\j Z 

believed that with the open ratio r^ = 0.0073 needed for zero interference at a small 
wing, the variation of 6 n at £ = 0 with y at a = 0 would be less than that shown 
in figure 6, since the variation is between 0.145 at x = 0 and 0 at x = 90°. Because 
of the image arrangement used to satisfy the assumed boundary conditions, which with 
y = 0, o' = 0 are no longer valid at the open boundaries, the values calculated at y = 0, 
a = 0 are believed to be slightly larger than would have been obtained if correct boundary 
conditions had been used. The value of 6 r x at 2 s o and a = 0 is zero for the 
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Interference factors, S R and 8 
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hovering condition x = 0 as it obviously must be. The fact that with ^ = 0 and x = 0, 
6 r is very close to zero for nonzero values of a is a little surprising, but no reason 
could be found to doubt the correctness of these calculated results. 

Although the skew angle x is treated as a parameter in the presentation of the 
calculated values of the interference factors, it may not be arbitrarily chosen, but depends 
theoretically through equation (C6) on the velocities V and Wq and on the angle a; and 
it increases toward 90° with increase in the stream velocity V more strongly than indi- 
cated by that equation because of the tendency of the wake to roll up into two trailing vor- 
tices like the vortex wake of a wing. Thus, reference 4 shows that for most operating 
conditions, the deflection of the wake from the horizontal is only about half that calculated 
by use of equation (C6). 

This behavior means that the uncertainty of the boundary conditions in the region of 
penetration of the open boundary by the helicopter wake can seriously affect the calculated 
interference in the neighborhood of the model only for a small range of stream velocities 
V near zero, where also x approaches 0. An indication of the variation with skew 
angle x °f the possible error arising from this uncertainty is seen in the contributions 
to the interference factors 6 r z and 6 R x at a = 0 due to subtraction of the part of 
the helicopter wake below the open lower boundary. Figure 7 shows that for the helicop- 
ter rotor mounted at the center of a test section with height-width ratio of 1.5, the effect 
of any such error becomes negligible at the rotor for x ~ a greater than about 45°. At 
the tail the effect becomes negligible for x ~ a somewhat greater than 45°. With height- 
width ratio less than 1.5, or otherwise with the rotor closer to the open lower boundary, 
the region for necessarily small error would be shifted to greater values of x ~ a - 



Figure 7.- Contribution to interference factors 6r z and 6r x due to subtracting part of helicopter wake below 

the open lower boundary, a = 0. 
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An investigation of the interference with a lifting jet in the test section closed on 
top, slotted on the sides, and open on the bottom was originally intended, but upon exami- 
nation, this problem did not seem to differ sufficiently from the helicopter interference 
problem to justify the additional labor of separate treatment. The essential difference 
is that the jet issues from the nozzle at some given angle a between its cross section 
and the tunnel test stream and with skew angle y equal to zero, but it is then swept 
downstream and thus the skew angle becomes a function of distance along the jet axis. 

The path of the jet could be estimated and segments of the jet could be represented by 
lengths of vortex cylinders at different skew angles to calculate the interference in a man- 
ner similar to the calculation for the helicopter. However, it is believed that the nature 
of the interference would be essentially the same as that of the helicopter and that the 
interference would in most of the operating range be not much different from that cal- 
culated for the helicopter at some average effective angle. This conclusion is partially 
substantiated by calculations given in reference 12. 

Other Configurations 

If the requirement for small variation of the interference with decreasing skew 
angle is relaxed (as it may be if the lift loading is small enough to prevent the wake being 
greatly deflected from the horizontal), the class of test-section configurations available 
for reducing the interference is considerably broadened. One which has been suggested, 
and which has the practical advantage of easy access to the test section, has completely 
open top and sides and solid bottom. For theoretically zero interference in a given test 
section, the model must be placed at a fixed distance from the solid wall. Appropriate 
distances above the lower boundary for zero interference z = 0) can be read from 
figure 8 for various ratios of height to width of the rectangular test section h/b. The 



Figure 8.- Interference factor 6r z for small rotor on vertical center line of test section 
of height h and width b with open top and sides and closed bottom as function of 
nondimensional distance g/h of rotor above bottom boundary. Wake horizontal. 


16 


values for h/b equal to 0.5, 0.66, and 1.0 have been derived from figures 14(a), 14(b), 
and 14(c), respectively, of reference 3 for the condition x ~ <* = 90°. For 1.5, the 
wing interference factor th# was computed as in equation (2) from the upwash velocity 
v for a lifting doublet computed as in appendix A except for the change in boundary con- 
ditions corresponding to replacement of the closed boundaries with open boundary and 
. replacement of the open bottom with closed boundary. The rotor interference factor may 
then be computed from the relation 

5 R,z* 46 w 

From figure 14 of reference 3 it is seen that the interference varies markedly with the 
skew angle x; thus, to maintain zero interference, the position of the model would have 
to be different for different values of x* Such variation of model position would some- 
what complicate model testing procedures. 

In wind-tunnel tests of lifting models for which the vortex wake is not deflected 
much (x - a *90°), the tunnel-boundary-interference problem is essentially that usually 
treated for lifting wings, and the interference velocity at the model can be reduced to zero 
by means of slots in all four walls as shown, for instance, in references 9 and 13. Such 
a configuration might be expected theoretically to produce more nearly uniform flow in 
the region surrounding the model than would the previously discussed test section with 
closed top, open bottom, and slotted side walls; but some concern might be felt regarding 
the danger of malfunction of the slots in the top wall because of the inflow of low- speed 
air in that region. Because of the panels in the bottom, the interference in this configura- 
tion should be more sensitive to skew angle for x - a near zero than would that in the 
configuration with completely open bottom. 

In an experimental investigation carried out by the NASA Langley Research Center 
using a high- lift wing, a test-section configuration with top closed and sides and bottom 
slotted seemed to produce more nearly uniform flow in the region of the model, particu- 
larly in the critical region of the tail, than did the previously suggested configuration 
with bottom actually or effectively (screened or many- slotted) open. This configuration 
was therefore investigated theoretically. Since a three-dimensional solution such as was 
previously obtained for the helicopter in the test section with open bottom would, with the 
slotted bottom, be mathematically complicated and time consuming to calculate, the two- 
dimensional method of reference 13 was applied. This method yields the interference 
velocity at the position of a wing (for x ~ a = 90°), but not that at the tail. Since this 
application of the theory of reference 13 is somewhat involved and since a part of the 
application to be made requires a small extension of that theory, the equations for cal- 
culation of the upwash interference factor Sw are developed in appendix D. 

Equation (D23) was used to compute the upwash interference factor 5^ at the cen- 
ter of a center- mounted wing represented by trailing vortices spaced 0.58 of the 
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test-section width in a rectangular test 
section with height 1.432 times the 
width, with closed top, with side walls 
and bottom having four equispaced slots 
each, and with (solid) panels cutting 
the axes through the test- section cen- ’ 
ter. The upwash interference factor 
is plotted in figure 9 against ratio of 
open area to total area of the slotted 
walls Tq. The interference at the 
center of the wing is indicated to be 
In the previously mentioned experi- 
ments approximately zero interference at the model was obtained with 10 percent open 
ratio in the sides and 15 percent open ratio distributed over six slots (close to an effec- 
tively open condition) in the bottom. The interference was still approximately zero when 
the bottom boundary was completely open. Thus, experimentally, approximately zero 
interference is obtained with slot openings more than four times as wide as indicated 
theoretically. This apparent divergence from theory may be at least partly due to the 
insensitivity of the theoretical interference to slot width in the range of width near to but 
greater than that width required for theoretically zero interference, but it also seems 
reasonable to suppose that with the narrow slots required for zero lift interference, vis- 
cous effects at the edges (and perhaps other departures from the assumptions, such as 
violation of the small-disturbance assumption) may result in boundary conditions approxi- 
mately corresponding in the theory to slot widths smaller than those actually being used. 

It is therefore suggested that not only should slot widths be experimentally designed to 
produce approximately zero interference for the types and sizes relative to the test sec- 
tion of models to be tested, but also consideration should be given to the possibility of 
Reynolds number effects on the slot performance. It is a fortunate circumstance that, 
because of the aforementioned insensitivity to width in the range of interest, slots that 
are wider than required, as they may be if they are designed in model-tunnel tests at 
relatively small Reynolds number, may still correspond to approximately zero interfer- 
ence. Moreover, this insensitivity to slot width favors compromise to obtain near-zero 
blockage interference, for which the zero-interference conditions require wider slots 
than required for zero lift interference. 

UPWASH INTERFERENCE AT THE TAIL 

It is evidently not very difficult, at least for a given model configuration, to design 
a test section for essentially zero upwash interference at the center of the wing, although 



Figure 9.- Upwash interference factor for center-mounted wing 
spanning 0.58 of width of rectangular test section having height 
1.432 times width. Top wall solid; side walls and bottom slotted 
with four equispaced slots each; trailing vortices horizontal. 


zero for an open ratio Tq of about 2.2 percent. 
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the open ratio Tq of the slotted sides may be different from that theoretically obtained. 

A much more difficult problem is to design a test section producing only small upwash 
interference throughout the region occupied by the test model. The upwash interference 
in the region of the tail is of particular importance. 

Because of the aforementioned mathematical complexity and time-consuming cal- 
* culations required, even with a large machine- computing system, the upwash interfer- 
ence at the tail has not been computed for the test section slotted on sides and bottom. 

» However, some enlightening qualitative considerations can be made. For this test sec- 
tion the upwash interference at the wing was computed as half the interference due to the 
interaction of the boundaries on the horizontal trailing vortex system at the far- 
downstream (Trefftz) plane. Since for the horizontal trailing- vortex system, the veloc- 
ities (zero far upstream) increase monotonically in the downstream direction, it follows 
that if the interference due to the trailing- vortex system without consideration of the 
boundary vortex is zero at the wing, it is zero also at the tail, provided the tail lies in 
the horizontal center plane of the test section. Some interference might result if the tail 
lies above or below the center plane, but this interference should be small in this type of 
test section because of the small vertical gradient of upwash velocity in the region near 
the horizontal center plane. Of greater concern is the effect of the downward deflection 
of the wake. 

In order to obtain a qualitative estimate of this effect, the two-dimensional theory 
of appendix D was used to calculate the upwash interference factor at the center of the 
test section for the wing displaced various amounts below the center. Such an estimate 
applies practically for relatively small wake deflection. For large wake deflection, for 
which the wake may even approach and penetrate the slots in the lower wall, a theory com- 
parable to the previously presented helicopter theory would be required. For this closed- 
top slotted- side-wall slotted-bottom configuration, the open ratio Tq was 2.2 percent, 
which is seen from figure 9 to correspond theoretically to zero interference at the center 
of the test section for the wing 
mounted at the center. From 
figure 10 it is seen that (as 
expected) movement of the wing 
(or rather of the trailing-vortex 
pair representing the wing) 
toward the slotted lower bound- 
ary results in a downwash inter- 
ference at the center of the test 
section. Downstream at the 
position of the tail of a center- 
mounted model, the downwash 
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3 Distance of wing below center^y/b 


Figure 10.- Variation of upwash interference factor at center as wing 
centered on vertical center line and spanning 0.58 of test section is 
moved below center. Height of test section 1.432 times width; top 
closed; side walls and bottom slotted with four equispaced slots each; 
open ratio r 0 = 2.2 percent; trailing vortices horizontal. 
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may be somewhat greater than (but less than twice) that indicated in figure 10. On the 
other hand, for the center- mounted model the vortex wake must start approximately in 
the horizontal center plane and in most cases should not be deflected much below the tail. 
It seems therefore that the downwash interference at the tail caused by deflection of the 
vortex wake in the test section with slotted sides and bottom should be small. With the 
completely open (that is, softer) lower boundary, the downwash at the tail due to wake 
deflection should be somewhat greater, and it is perhaps for this reason that the slotted- 
bottom configuration seemed experimentally to produce smaller interference at the tail. 

In contrast to the interference at the wing, that at the tail requires consideration 
of the interaction of test- section boundaries on the bound vortex as well as on the trailing 
vortices; that is, the estimation of interference at the tail is a three-dimensional problem 
For the closed-top open-bottom configuration, a qualitative consideration is easy. The 
top- and bottom -boundary conditions are satisfied by image vortices. The first three 
images above and below the original bound vortex are indicated in figure 11, where the 
sense of the vortex is in each case indicated by an arrow. The vortex images 1 combine 
as do also vortex images 3 to produce zero upwash on the horizontal center plane. A 
small upwash is produced by images 2, and even this upwash is somewhat reduced by the 
series of images at distances greater than three test-section heights above and below the 
original bound vortex. Note that this upwash tends to cancel with the downwash caused 


Image 3 “ “O 


(mage 2 — 


image I-— Qj 

^ — Closed boundary 

Original 

bound vortex " ^ 

rjItiL" x 
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— Lower boundary 

Image 2 *— 


Image 3 — (^y 

® Induced velocities on axis due to images 1 

® Induced velocities on axis due to images 2 

@ Induced velocities on axis due to images 3 


Figure 11.- Illustration of upwash interference of closed top and open bottom 
boundaries in presence of the bound vortex. 
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by deflection of the trailing-vortex wake. For a slotted rather than open lower boundary, 
the resulting upwash at the tail due to interference on the complete vortex system would 
be slightly increased. 

To the extent that it is allowable to treat the bound vortex separately from the 
remainder of the vortex system, solid side walls produce no interaction, inasmuch as 
* vertical planes normal to the vortex axis contain streamlines of the flow. At slots, how- 
ever, the constant-pressure boundary condition requires a flow opposing the vorticity, 

. and for the lifting vortex this flow is upward at the tail position. This effect should be 
small as long as the span of the lifting device is small relative to the test- section width. 
From this qualitative theoretical consideration, it is seen that for a zero- lift-interference 
(at the wing) test section with closed top, slotted side walls, and slotted or open bottom, 
and with the lifting wing mounted at the center, 

(a) The interference at the tail due to interaction of the boundaries with a horizontal 
trailing-vortex system is zero provided the tail is located near the horizontal center line, 

(b) Passage of the vortex wake horizontally somewhat below the tail produces down- 
wash interference at the tail, and 

(c) Interaction of the test- section boundaries with the bound vortex produces upwash 
interference at the tail. 

These conclusions are, of course, conditioned by the assumption of potential flow at the 
slots. 

In order to obtain further information concerning the contribution of slotted side 
walls to the upwash interference at the tail, equation (C61) was used to compute upwash 
interference factors for a lifting element centered between two vertical-plane infinite 
parallel slotted walls spaced a distance b apart. The development of equation (C61) 
implies application to a small, lightly loaded wing (horizontal wake) for which the inter- 
ference on the bound part of a possible vortex-system representation is included. The 
stream flow (velocity V) is assumed to extend infinitely far above and below the wing. 
Examination of the right-hand side of this equation shows that at the position of the wing 

= o), the interference is given entirely by the second term. On the horizontal cen- 
ter line far downstream ^ ^ = 0^, on the other hand, the first term makes a contri- 

bution equal to the second; thus, in agreement with a well-known theorem, the upwash 
interference at the wing is exactly half that at the far- downstream (Trefftz) plane. 

The value of the restriction constant l (see eq. (B2)) required for zero upwash 
interference at the wing was found to be 0.592b. This value is about the same as that 
needed for zero interference in the test section with closed top and open bottom and with 
height 1.432 times its width. Some consideration shows the nature of the slotted-wall 
effect. The flow through the slots far from the wing permits a general downflow to 
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Distance above horizontal center line,z/b 


Figure 12.- Variation with distance above horizontal center plane of upwash interference factor for lifting element 
centered between infinite vertical parallel slotted walls spaced distance b apart, l/b = 0.592. Wake horizontal. 







develop in the whole field above and below the wing, an effect similar to the downflow 
known to develop in a test section having closed sides and open top and bottom. The open 
ratio Tq required for zero upwash interference must be of such a value that at the wing 
this downflow is exactly counterbalanced by the upflow induced by the vorticity developed 
in the presence of the lifting wing by the portion of the side walls opposite the wing. Com- 
puted upwash interference factors in the vertical center plane above and downstream from 
the wing are shown in figure 12. These results are in agreement with the qualitative 
analysis here presented. The interference velocities, which at 2 = 0 are zero at the 
wing, show the expected downwash above the wing. Since equation (C61) is even in z/b, 
exactly the same variation occurs below the wing. Superposed on this variation above 
and below the wing is a much smaller streamwise variation (the interference corre- 
sponding to the bound part of the vortex system) which, as seen from the first term on the 
right of equation (C61), is odd and therefore for upstream (negative) values of x/b is 
the negative of that indicated in figure 12. 

The downstream variation of the downwash interference factor on the horizontal 
center line is shown in figure 13. The previously deduced upwash in the region likely to 
be occupied by the tail is evident. The interference velocity on the center line reduces 
to zero for x/b approaching °° (for zero interference at the wing) but at ^ = 0.6 
it is still increasing. 

The general downwash interference induced with the parallel infinite walls shows 
the importance of the top and bottom walls (perhaps also of closed-test-section entrance 
and exit) and suggests that the uniformity of flow in the region of the wing is favored by a 
moderate ratio of test- section height to width (for example, less than 2). 

Note that the upwash interference factor shown in figures 12 and 13 is a helicopter 

rotor interference factor based on the square of the width b and therefore must be 
1 432 

multiplied by — — before comparison with the upwash interference factor for the 
wing in the test section with height-width ratio of 1.432. 
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For the test section with closed top wall, slotted side walls, and open lower bound- 
ary, the interference at the tail can be calculated from the theory of appendix C. The 
upwash interference factor 6^ z along the tunnel center line downstream from a lightly 
loaded lifting element (small wing, a = 0, % = 90°) mounted at the center of such a test 

section with height-width ratio ^ = 1.5 was calculated from equation (C65). The calcu- 

b 

lated values for r^ = 0.0455, as in figure 6, and for Tq = 0.0073, corresponding to zero 
upwash interference at the lifting element, are presented in figure 14. For r^ = 0.0455, 
a downwash interference occurs at the wing and the downwash increases monotonically 
downstream. For Tq = 0.0073, the upwash interference is zero at the wing; downstream, 
it is due entirely to the interference on the bound vortex and, as might be expected and as 
indicated by the qualitative theory previously presented, after first increasing it decreases 
toward zero at large distances downstream. 

The values of 6^ z on the upper curve of figure 14, which corresponds to x~ 90°, 

may be compared at ^ * 0, 0.1, 0.3, and 0.5 with the values for other skew angles x 
at a= 0 in figure 6. The interference is seen to vary from upwash to downwash as x 
is increased from zero toward 90°, and the variation increases with distance jat least 

to ^ = 0.5^ downstream from the lifting element. At ^ = 0.5 the downwash reaches 
values near x = 70° significantly greater than that indicated in figure 14. This varia- 
tion cannot be explained as being caused by displacement of the wake as in figure 10, but 
is evidently a function of the wake angle itself as well as of the position of the wake rela- 
tive to the field point of interest and to the lower boundary of the test section. 



Distance downstream from wing,x/h 

Figure 14.- Upwash interference factor along tunnel center line downstream 
from small wing mounted at center of test section with closed top wall, 
slotted side walls, and open lower boundary, h/b = 1.5; wake horizontal. 


24 



CONCLUDING REMARKS 


By use of a test section having top wall solid, side walls slotted with four equispaced 
slots each, and bottom boundary open, it is found possible to reduce to zero the theoretical 
wind-tunnel-boundary upwash interference at a small, lightly loaded (trailing horizontal 
' wake) wing mounted at the center of the test section. For a high-lift-coefficient model, 
for which the wake is not horizontal, the theoretical upwash interference in this test sec- 
tion does not vary greatly with variation of wake deflection. The slot width required for 
zero lift interference at a wing is found to be very small, but because of boundary -layer 
and viscous effects, as well as other uncertainties in practical application of the theory, 
the theory is regarded as unreliable for determining the slot width, and an experimental 
determination is therefore suggested. It is a fortunate circumstance that the lift inter- 
ference at a small wing is relatively insensitive to slot width for widths near and some- 
what greater than that required for zero interference. Aside from easing the problem 
of slot-width selection, this behavior permits compromise in favor of reducing solid 
blockage interference, which is known to require wider slots than are needed to produce 
zero lift interference. 

For a test section of the type herein suggested, the theoretical interference upwash 
on the tunnel center line downstream from the lifting element is slightly positive, pro- 
vided the slot width is such as to produce zero upwash interference at the wing and pro- 
vided the vortex wake from the wing lies in the horizontal center plane. However, there 
is reason to believe that even with this ideal arrangement, a downwash interference 
exists both above and below the horizontal center plane so that if the tail is appreciably 
off this horizontal plane, it may experience a downwash rather than an upwash. This 
effect, along with the effect of movement of the vortex wake toward the open lower bound- 
ary and the effect of wake angle itself, leads to an appreciable downwash at the tail if the 
wake is even moderately deflected, for example, 10°. Moreover, if the slot width is 
greater than that needed for zero lift interference at the wing, the corresponding down- 
wash produced downstream from the wing is greater than that at the wing. In practical 
application, therefore, the interference at the tail of a winged model mounted at the center 
of such a test section is very likely to be of the nature of a downwash. On the other hand, 
for high- lift- coefficient models, the interference on the tunnel center line downstream 
from the position of the model tends toward an upwash as the vortex wake approaches the 
downward vertical direction. 

Zero interference at a small wing with horizontal vortex wake can be obtained with 
other test- section configurations, but to prevent large variations with skew angle while 
maintaining small interference at the model, an open or effectively open (many-slotted or 
screened) lower boundary seems necessary. However, the use of such a boundary intro- 
duces additional uncertainties in boundary conditions and therefore further degrades the 
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reliability of the interference calculations relative to such calculations for completely 
closed boundaries. A many- slotted or screened, rather than completely open, lower 
boundary may be desirable to prevent oscillations of the test-section flow. If the require- 
ment for small interference variation with skew angle can be relaxed, the zero interfer- 
ence on a wing with horizontal wake can be obtained in a test section having all four walls 
slotted with coarse slots. The flow in the region of the wing might then be more nearly 
uniform than that for the test section herein suggested, and the effect of moderate down- 
ward deflection of the wake might also be less. With high-lift devices a slotted upper 
wall might behave essentially as an open boundary because of the effectively thick bound- 
ary layer produced by inflow through the slots, but there is reason to believe that this 
behavior would not affect the interference much provided the bottom wall were suffi- 
ciently effectively rigid to carry the reaction to the lift. In contrast, with an open lower 
boundary a too great lift loading on the test- section flow might result in flow separation 
from even a solid upper wall, and the resulting interference would be a downwash as in 
a test section with completely open top and bottom walls. 

Because the practical boundary conditions at mixed (open and closed) boundaries 
are not well determined, particularly with high lift and narrow slots, the theory herein 
presented is regarded as unreliable for predicting numerical values for lift interference 
in given test- section configurations. It should be useful in showing the nature of the 
interference and the possibility of reducing it, and it should serve as a guide to design 
and experimentation. The slot widths required for zero interference should be deter- 
mined by experiment, although the experimental investigation of wind-tunnel-boundary 
interference is difficult, largely because the quantities of interest are usually relatively 
small differences. In a zero- interference (at the wing) test section the distribution of 
upwash velocities and their variation with skew angle might well be approximated by the 
theory even though the practical slot widths are different from those theoretically deter- 
mined. It is believed that further theoretical investigation might lead to test-section con- 
figurations having more nearly uniform flow in the vicinity of the model, including the 
position of the tail, than exists for those so far investigated. With considerable additional 
labor, the effects of curvature of the wake (as for downward-directed jets) might also be 
taken into account. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 1, 1968, 

721-01-00-20-23. 


26 



APPENDIX A 


DERIVATION OF EQUATION FOR CALCULATING UPWASH INTERFERENCE 
VELOCITY ON SMALL LIFTING MODEL IN TEST SECTION WITH 
CLOSED TOP AND SIDE WALLS AND OPEN LOWER BOUNDARY 


The model is located on the vertical center line of a rectangular test section of 
height h and width b with closed top and side walls and open lower boundary. The 
model is located at height g above the lower (open) boundary. The upwash interference 
velocity v is that induced by the test- section boundaries on a vortex doublet corre- 
sponding to the lift of the model and trailing horizontally downstream from the model. 

The induced upwash velocity is first obtained in a plane normal to the vortex doublet line 
far downstream (the Trefftz plane), where the flow is essentially two-dimensional. The 
induced upwash velocity at the model is then just half that at this far-downstream plane. 

The symbols peculiar to this appendix are as follows: 


n,k any integers 

yj,Zj Cartesian coordinates of a vortex doublet 

<f> + velocity potential of line doublet corresponding to positive lift 

<p_ velocity potential of line doublet corresponding to negative lift 


The arrangement of vortex doublet images satisfying the boundary conditions is 
indicated in figure 1. With the origin of y,z coordinates taken at the position of the 
original lifting doublet, vertical image rows occur at y * kb, k being all integers from 
- 00 to °°. Horizontal rows of positive (lifting) vortex doublets occur at z = (2n)(2h) 
and at z = (2n)(2h) - 2g, n being all integers from -°° to 00 . Horizontal rows of 
negative vortex doublets are located at z = (2n + l)(2h) and at z = (2n + l)(2h) - 2g, 
n being all integers from -°° to °°. 


The velocity potential of a lifting-line vortex doublet of strength n 
is (see ref. 14) 

, M z 
+ 2ir y2 + z 2 

The potential for the line doublet of opposite sense is 


located at (0,0) 
(Al) 


= 


-M z 

2 y2 + z 2 


(A2) 
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If the vortex doublet is located at y^Zj, the potential is 

JL 

2l7 


0 = JL Llfl 

5- - -2 . „ ,2 


(y - yi) + ( z - z i) 


or 




Z - Z- 


(y - yi)‘ + (z - z i)‘ 

On the Z-axis, y * 0 and these potentials give 

z - Z-, 


0 =-^ 

+ 2 77 


( Z - Zj)2 +yi 2 


(A3) 


(A4) 


(y = 0) (A5) 


0 «=lii 

2tt 


Z ” Zi 


( z • z i ) 2 + n 2 


For z * zj the sum for any horizontal row is 


^(row)+ it Z 


Z - Z, 


or 




k=-°° 


00 

1 


( z - 


2 + k 2 b 2 


Z - Zi 


(y = 0) (A6) 


(y = 0) (A7) 


(y = o) (A8) 


(row)- " 277 ^ ( z _ Zl )2 + k 2 b 2 
As indicated by equation (I) of reference 8, the sum of the infinite series in equation (A7) 


is 


h , = JL ooth 

v (row)+ 2b 


b 


(y = 0) (A9) 


and that of the infinite series in equation (A8) is 


^(row)-=ib COth ^ b'"" 1 ' 

The corresponding upwash velocities at y = 0 are 


(y = 0) (A10) 


V/ * - JL 

(row)+ 2b 


- i csch 2 ± - z 0 
b b 


= - csch 2 


ff ( z ~ z l) 


(y = 0) (All) 
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V (row)- 


= csch 
2b 2 




(y = 0 ) (A 12 ) 


In order to find the interference, the potential of the doublet located at 0,0 must be 
removed from that of the row at z^ = 0 . Since the potential is required on the Z-axis in 
the neighborhood of z = 0, the potential may be expressed as 


0 (row) = coth = 


1 1 Arz\ _ 1 firz\ 

7rz/b 3 \b ) 45 \ b / 



(A 13 ) 


(y = 0; zj = 0) 

Hence the interference upwash at 0,0 due to the horizontal row at z^ = 0 is 7rju/6b 2 . 

At the position of the model the upwash is just half that at the lar downstream Trefftz 
plane so far considered. The total interference upwash at the position of the model where 
z approaches 0 is therefore 


v = JDL. _ UL csch 2 M - JUL S csch 2 ^ + csch 2 47r - nh ~ 2ff S 
12b 2 4b 2 b 4b 2 n J - J ^ \ k b < 

n^0 


oo 

y 

4 b 2 ^ 


n=-oo 


csch 2 — (^-i~- 1 ) h + csch 2 fo( 2n + *) h 
b b 




With some rearrangement this equation can be rewritten as 


v = ^2 ■ 4 ■ csch2 ^-|( 2 csch2 nf + csch2 + csch2 -- h b " ?g ) 


°° r- 


n=l 


2 csch 2 27r < 2n - ^ + csch 2 MglLl Dh - 27 rg + cgch 2 27 r( 2 n - l)h + 2 ff gl v (A14) 
b b b ' 


This form is easily seen to be an alternating convergent series of monotonically 
decreasing terms. The error in cutting off the summation at any term is therefore less 
than the value of that term. The convergence is so rapid that only the first few terms 
need to be considered. 
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DEVELOPMENT OF EQUATION FOR CALCULATING UPWASH INTERFERENCE 
IN TEST SECTION WITH CLOSED TOP, SLOTTED SIDES, 

AND OPEN LOWER BOUNDARY 

The configuration is the same as that of appendix A except that the side walls are 
now slotted; therefore, the boundary condition at the slotted walls must be satisfied 
directly rather than by means of images. The symbols for this appendix are as follows: 

A(co) arbitrary function of co determined to satisfy slotted boundary condition 

b width of test section 

d distance between centers of two adjacent slots 

g distance of model from open lower boundary 

h height of test section 

l restriction constant defined in equation (B2) 

n any integer 

q variable of integration, cob 

Tq ratio of slot width to distance d 

V tunnel stream velocity 

v upwash velocity at the model 

y,z rectangular Cartesian coordinates 

yj,Zj coordinates of a line vortex doublet 

A6i contribution to upwash interference factor at 0,0 due to a line vortex image 
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A6 S contribution to upwash interference factor at 0,0 due to interference of 

slotted side walls on a line vortex doublet 

6w upwash interference factor 

/i strength of vortex doublet 

4> velocity potential 

<t>£ velocity potential of vortex doublet 

<p s velocity potential induced by slotted side walls in presence of vortex doublet 

transform of (f> on z - zj 

O s transform of (p s on z - 

co variable of transformation on z - z^ 

For simplicity the slotted- side-wall boundary is replaced by a homogeneous bound- 
ary, the boundary condition on the velocity potential (p being given by (see ref. 9) 

4> ± l = 0 (Bl) 

9y 

where the positive sign applies at the boundary y = ~ and the negative sign applies at the 
b ^ 

boundary y = - -. Equation (Bl) is derived under the small -disturbance assumption for 
slots in the direction of the tunnel stream V; thus, the velocities at the boundaries nor- 
mal to the slot direction are small. The restriction constant l is given by equation (3) 
of reference 9 as 

H 7rr n 

l = f log e esc (B2) 

where d is the distance between the centers of two adjacent slots and Tq is the ratio 
of slot width to the distance d, or for a uniformly shifted wall simply the proportion of 
the wall that is open. The restriction constant has the dimension of a distance and 
approaches zero for an open tunnel and infinity for a closed tunnel. The homogeneous 
boundary approximation is satisfactory for determining the interference at distances 
from the wall that are large relative to the distance between slots. (See, for example, 
ref. 9.) 
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The boundary condition (Bl) is satisfied by adding to the potential 4 > ^ of a vortex 
doublet located at 0,zi in the far-downstream plane (Trefftz plane) normal to the test 
stream a potential <i> s such that + 0s satisfies equation (Bl). Of course, 4> s as 
well as $d must satisfy the Laplace equation 

+ = 0 (B3) 

8y “ 3z“ 

From equations (A3) and (A4) the potential of a positive (lifting) line vortex doublet loca- 
ted at 0,zj is 

A. = JL Llfl (B4) 

d>+ 2 tt y2 + _ Z jj2 

and the potential of a negative line vortex doublet located at 0,zj is 





(B5) 


In order to obtain the additional potential <j) s due to the influence of the slotted 
side walls, the exponential Fourier integral transform (see ref. 15) on z - zj with 
variable of transformation to is taken. Transformation of equation (Bl) then gives 


ft(y,co) 


± 1 8 »( y >^ ) = o 

a y 


(B6) 


The Laplace equation (eq. (B3)) transforms to 

8y2 


(B7) 


A solution of equation (B7) intended to correspond to </> s and having the required sym- 
metry in y is 

S2s(y,w) = A(co)cosh(coy) (B8) 


where A(co) is an arbitrary function of to. The potentials (eqs. (B4) and (B5)) of the 
line doublets transform to 




-MM 


(B9) 


This transformation is easily derived from formula (15), p. 65 of reference 15. 

Substitution of for in the boundary condition (eq. (B6)) at the right- 

hand wall y = x gives 
& 
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cosh^o* + loo sinh (co ^ 


-WfSr-^ |H 


AM 


The same equation holds at the left-hand boundary. Consequently, 


AM = ± 


Sr te )‘ 


■|"i 


cosh + loo sinh kj 


and the transform S2 s of the interference potential </> s due to the reaction of the slotted 
side walls on the line doublet potential 0^ ± is 


= ± ^ 


iW 


cosh^o) + loo sinhfo) ^ 


cosh(coy) 


(BIO) 


where the positive sign applies for the original lifting doublet or its positive image and 
the negative sign applies for its negative image. 

By inversion, 


P oor 


‘ * u 


O 


(0 


-rM 


Me 2 


cosh^cu + loo sinh^co | j 


cosh(cuy) 


iwfz-zi ) 

5 V 1/ l 


dco 


The corresponding contribution to the upwash velocity at y = z = 0 is 


(Bll) 


d(p 


s,± 


8z 


y=o 

z=0 


= *f 

4tt 


CO 

Icol 


- loo e 


2\ -jH 


cosh^c 0 |j + loo sinh(a> ^ 


-1G0Z1 

e do? 


27 T 


-a& 


(<■» - iJk "2 


0 


cosh 


H 


+ ico sinh (u> ^ 


cos 


(cozi) 


dco 


(B12) 
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Jib 

Since, as given by equation (2), the upwash interference factor is 6^ « — v, where 

v is the upwash velocity at the model (half that in the Trefftz plane), the contribution to 
the upwash at 0,0 is 


A6 


s,± 


x hb 


(a. - *<, 2 )e“"2 

j 

0 

cosh^cu + icu sinh^cu ^ 


rs<x 

> _ cub 

= T i_l 1 


|cub - ^ cu^b^je ^ 

8t t b 

0 

cosh/^W ^ cub sinhf^ 
\ 2 / b \ 2 


cosfcuz^dcu 



and with cub * q, 


A5 


s,± ~ 


=F 


_1_ h 
8irb 


f>oO 


JO 



cosh x + r Q sinh ^ 
2 b 2 


cos 



dq 


(B13) 


From equations (B4) and (B5) the upwash velocity at 0,0 due to a doublet image at 
0,zj is found to be 

8 »d, ± 


9z 


■ = T 


27TZJ 2 


and therefore the contribution to the interference factor at the model is 

A 6, 


= T 1 b 1 
%± 87rh/„.\2 


\h / 


(B14) 


As for the configuration with solid top and side walls and open bottom boundary, positive 
images in the top and bottom boundaries occur at zj = (2n)(2h) and at zj = (2n)(2h) - 2g 
for all positive, negative, and zero values of n, the original lifting doublet being at zj = 0 
Negative images occur at z\ = (2n + l)(2h) and at z\ = (2n + l)(2h) - 2g where n is 
all integers from to °° and g is the distance of the model from the open lower 

boundary. It follows therefore from equations (B13) and (B14), since the velocity due to 
the original lifting doublet is not to be included, that the upwash interference factor at the 
center of the tunnel is 
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(B15) 


Convergence should be rapid; thus, only a few positive and negative values of n are 
needed. 
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DERIVATION OF EQUATIONS FOR CALCULATING BOUNDARY- INTERFERENCE 
FACTORS FOR A HELICOPTER ROTOR IN A RECTANGULAR TEST SECTION 
WITH CLOSED TOP, SLOTTED SIDE WALLS, AND OPEN LOWER BOUNDARY 

A lifting helicopter rotor is assumed to be mounted at the center of a rectangular 
wind-tunnel test section with closed top wall, slotted side walls, and open lower boundary. 
Equations for calculating the streamwise interference factor x and the upwash inter- 
ference factor 6j^ z induced by the test-section boundaries are to be derived. Symbols 
applicable particularly to this derivation are as follows: 

A area of rotor disk 

A vector area of rotor disk 

B(w,g) arbitrary function of w and g 

B s (w,g),Bg(w',g') values of B(w,g) appropriate to dG ns and dG^ s , respectively 
b width of test section 


c 

cross-sectional area of test section, 

hb 


G 

exponential Fourier transform of 

ft 

on z 


G n 

exponential Fourier transform of 

ft n 

on z 

- 2nh + s cos(x - Qf) 

Gn 

exponential Fourier transform of 

«n 

on z 

- (2n - l)h - s'cos(x - ot) 


G n ,s>G^, s values to be added to G n and G^ to satisfy slotted boundary condition 

g variable of transformation on z - 2nh + s cos(x - a?) 

g' variable of transformation on z - (2n - l)h - s'cos(x - <*) 

h height of test section 

i=fl 
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I 1’^2’ I 3» I 4 integrals defined by equations (C48) 

[I W1 ] , |lw 2 j contributions to potential of interference velocities arising from 

n = o n=0 subtraction of vortex cylinder extending below lower open 

boundary 


fj,k unit vectors in direction of X-, Y-, and Z-axes, respectively 

Kq,Kj modified Bessel functions of second kind 


restriction constant 


M}, M 2 parameters defined by equations (C 37) and (C42) 

n any integer 

P field point 

p = wh or cob 


q = gh or gb 


r distance from element of surface to field point P 

? vector distance from element of surface to field point P 

s distance from rotor disk to vortex ring 

T thrust 


u,v interference velocities in x- and z-direction, respectively 

V tunnel stream velocity 

Vr magnitude of resultant of V and Wq regarded as vectors 

Wq induced velocity at rotor disk 

X,Y,Z coordinate axes 
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x,y,z 

z o> z 6 

a 

r 

( A6 R,x)ns 


coordinates 

variables defined by equations (C38) and (C44) 
angle of attack of rotor-tip-path plane 
circulation, positive for positive lift 

slotted- side-wall contribution to streamwise interference factor at (x,0,0) 
for a vortex cylinder cell nominally located at (0,0, 2nh) 


( A ®R,z) ns 


slotted- side-wall contribution to upwash interference factor at (x,0,0) for 
a vortex cylinder cell nominally located at (0,0, 2nh) 


e small positive number 

6 angle between A and ? 

p density of test medium 

<p velocity potential 


velocity potential of vortex 


cylinder cell nominally lying at (0,0,0) 


<Pn velocity potential of part of vortex cylinder cell image corresponding to 

part of original vortex wake between rotor and open boundary 


$0,n velocity potential of vortex cylinder cell nominally located at (0,0, 2nh) 

4> w velocity potential of part of rotor wake outside position of open lower boundary 

<f >0 velocity potential of part of rotor wake outside position of open lower 

boundary plus image in that boundary 

< contributions to <p g given by equations (C41), (C43), (C45), and (C46), 
respectively 

X skew angle, positive counterclockwise 

a variable of transformation on x - s sin(x - a) 
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Cl exponential Fourier transform of <p on x 

£2 n exponential Fourier transform of <fi n on x - s sin(y - a) 

Primes on symbols indicate values applying to image of vortex cylinder reflected 
across open boundary. Subscript s indicates a slotted-wall contribution. 

As in references 2, 10, and 11, the disturbance caused by the helicopter rotor is 
represented by a straight elliptic cylindrical vortex sheet extending from the rotor. 

The vortex sheet is made up of a continuous constant- strength density distribution of cir- 
cular vortex rings lying in planes parallel to the rotor -tip-path plane. (See fig. 5.) In 
figure 5 right-handed axes X,Y,Z are chosen, but for simplicity the Y-axis is not shown; 
this omission causes no difficulty, because the rotor disk and the flow field are symmetric 
about the X,Z plane. The wind-tunnel stream velocity V lies in the direction of the 
positive X-axis. The rotor disk (tip-path plane) makes an angle a with this direction 
and is positive in the sense indicated in figure 5. The angle between the normal to the 
rotor disk and the axis of the cylindrical vortex sheet is designated y. 

The element of potential dc/> at the field point P due to the vortex ring of strength 
— ds from reference 16, page 212, is, in the notation of this paper, 

(C1> 

where r is the distance from an element of surface enclosed by the vortex ring to the 
point P, and s is distance from the rotor disk to the vortex ring. In the calculation 
of the boundary interference the points P of interest are considered to be far from the 
vortex rings contributing to the interference so that the relative variation of r over the 
surface A enclosed by any vortex ring is small. Since also the angle 6 between the 
vector area A and the vector distance ? is almost constant for the integration, equa- 
tion (Cl) is well approximated by 

a dr 

d <p = -^r cos 6 ds (C2) 

4nr 2 ds 

If A and 6 have the directions indicated in figure 5, f lies along A if 6 is zero. 
With r positive as indicated in figure 5, the velocity through the vortex ring is opposite 
to r in direction and corresponds to that of a lifting rotor. The partial derivative of 
d(p with respect to r gives a negative velocity (for 6 = 0) as should exist on the axis 
of the vortex ring for the lifting rotor. Thus, the sign in equation (C2) is correct. 

In order to express r and cos 6 in terms of x, y, z, and s, note that from fig- 
ure 5, with i,j ,k unit vectors in the directions of the X-, Y-, and Z-axes, respectively, 
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and 


so that 


and 


It follows that 


A = i A sin a + kA cos a 
r » i|x - s sin(x - a >] + Jy + k[z + s cos(x - a)j 

r 2 = £x - s sin(x - a)J 2 + y 2 + Jz + s cos(x - o?)J^ 

A . y [x - s sin(x - a)]sin a [z + s cos(x - Qf)]cos a 
~ Ar r + r 

dr 

^ ds^ ds [x - s sin(x - a)] sin a + [z + s cos(x - a;)] cos a 

d d> = - =! 

477 r 3 


cos 6 


(C3) 


In reference 10 it is shown that the density dr/ds of vortex strength along the 
cylindrical vortex sheet is numerically equal to the final induced velocity along the axis 
of the cylinder. It is also shown that this final induced velocity is numerically equal to 
the final induced velocity normal to the rotor disk. Moreover, the normal induced veloc- 
ity at the disk is half that in the (final) far-distant wake. Thus, dr/ds may be replaced 
with 2wq, where Wq is the normal induced velocity at the rotor disk, positive for posi- 
tive lift: 

g=2w„ <C4) 

Equation (C3) may therefore be written as 

Awg [x - s sin(x - a)] sin a + [z + s cos(x - a)"]cos a 
d<p = — — ^ ^ r-= ds (C5) 


Although the skew angle x and the velocity wg induced at the rotor disk are the sig- 
nificant parameters of the helicopter wake, they cannot be specified independently of the 
wind-tunnel stream velocity V. According to equation ( 8 ) of reference 10, the param- 
eters are related, in the notation of this paper, by 


V_ 

w 0 


2 tan ^ 
cos(x - a) 


(C 6 ) 


(Note that the angle a is here the negative of the angle designated a in ref. 10.) For 
skew angles other than zero, the induced velocitiy varies over the rotor disk, but wg 
may be regarded either as the induced velocity at the center or as the average of the 
induced velocities from front to rear over the diameter, since according to reference 10 
these values are equal. The average induced velocity wg may be used in the Glauert 
relation to compute the thrust T: 
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T = 2pAV R w Q 

where V R is the magnitude of the resultant of V and Wq regarded as vectors. 

In free air the vortex cylinder of figure 5 extends far below the lifting rotor; there- 
fore, the velocity potential is found by integrating the expression (C5) over s from 0 
to 00 . In the wind tunnel with open lower boundary, the rotor wake is assumed to pass 
out of the test stream and to have no further effect on the model. The possible conse- 
quences of this assumption have already been pointed out. Under this assumption part 
of the boundary interference is therefore represented by subtracting the potential (or the 
corresponding induced velocities) due to the part of the vortex cylinder extending below 
the open boundary as indicated by the dashed lines in figure 15. 

The zero-pressure-disturbance condition at the open lower boundary, which corre- 
sponds to zero induced velocities in the direction of the tunnel stream at that boundary, 
is satisfied by a mirror reflection of the vortex cylinder across the boundary without 
change in sense of the vorticity. The original disturbance and its reflection across the 
open boundary make up a vortex cylinder cell. The solid-boundary condition is satisfied 
by a similar reflection, but with reversal of sense of the vorticity. The complete image 
system satisfying the boundary conditions at the upper and lower boundaries in a test 
section with solid upper boundary and open lower boundary and with the rotor at the cen- 
ter of the test section consists of a vertical row of vortex cylinder cells spaced at twice 
the test-section height h along the Z-axis as shown in figure 15. The image of the 
original disturbance reflected across that boundary is of the same sense as that in the 
subtracted part of the original vortex cylinder and may therefore be expected to provide 
some compensation for the loss. However, the upstream displacement of vorticity may 
result in distortion of the flow downstream from the rotor. 


For a vortex ring lying on the image of the original disturbance vortex cylinder 
reflected across the open boundary, the element of potential (eq. (C5)) must be replaced 
with 


d <p = 


where 


Awq J 

r -l 

x - s’sin(x - a)] 

sin a + [ 

z + h - s’cos(x - 4 

cos a 

277 1 

r'2 

L 

[x - s'sin(x - a) 

) ^ 2 

] +y^ + 

3 

2 

jjz + h - s'cos(x - a)} 


ds' (C7) 


s' is the distance along the axis of the image in the open boundary starting from the point 
(0,0, -h), and h is the test-section height. 


The potential due to the original disturbance vortex cylinder and its image in the 
open boundary is obtained as the sum of the integrals of expressions (C5) and (C7), where 
s runs from 0 to its value at the open boundary and s' runs from 0 to its value at the 
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open boundary. This potential, designated 0 q, applies to the vortex cylinder cell nomi- 
nally located at (0,0,0). The potential of the vortex cylinder cell nominally loca- 

ted at (0,0, 2nh) is obtained by obviously similar integrations of 


n Awo ( [x - s sin(x - of)] sin or + [z - 2nh + s cos( x - a ) cos cn 
d^ n =(-l) ~ — I ~ — ' ds 


2 % 


S 7 s sin(x - a)J^ + y^ + - 2nh + s cos(x - a)]^j 


and 


, n Awq / -|x - s’sin(x - afljsin a + jz - (2n - l)h - s'cos(x - «)j cos a\ 


& 4 > n = (-1) 


2tt 


|^jx - s’sin(x - + y^ + jz - (2n - l)h - s'cos(x - 


(C8) 


ds’ (C9) 


where n takes on all positive and negative (as well as zero) integer values in the sum- 
mation of <p n over n to obtain the interference potential. However, because of the 
alternating signs only a few terms of the infinite summation on n should be required 
for practical convergence. To obtain the interference, the potential of the original dis- 
turbance vortex cylinder must be subtracted out of the summation. Thus, the potential 
of the whole semi-infinite vortex cylinder extending from the rotor must be subtracted. 
This potential is obtained by integrating expression (C5) on s from 0 to 00 . 

The boundary condition at the slotted side walls must now be satisfied. It is suf- 
ficient to satisfy the side-wall boundary condition separately for every original or image 
differential element d0 n and d <p' n of the disturbance potential. The side-wall bound- 
ary condition is as in appendix B (eq. (Bl)) : 

(CIO) 

b 

where the positive sign applies at the boundary y = k- and the negative sign applies at 
the boundary y = - ^ and where l is the restriction constant defined in appendix B 
(eq. (B2)). 

To facilitate the solution, exponential Fourier transforms of d(p n on 
x - s sin(x - a) with variable of transformation co and of d$^ on x - s'sin(x - «) 
with variable of transformation co T are taken. Then the following equations, corre- 
sponding to expressions (C8) and (C9), are given by formulas 7, page 11, and 27, 
page 66, of reference 15: 
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(-l) n Aw 0 ds / 

-i(sin oi)u)Ko 


d« n = 


{u\jy 2 + [: 


z - 2nh + s cos(x 




+ [jz - 2nh + s cos(x - a!)J(cos a) 


a>Ki| 


[y 2 

+ | 

\z - 2nh + s cos(x - 

«)] 


1 

F 

i | 

1 + 

[- 

- 2nh + s cos(x - a)] 

i* / 


+ [z - 


2nh + s cos(x 


- <f > o) 


(Cll) 


d^ = 


(-l) n Aw n ds' / 

’ U i(sin a)WKo 


jw’^y 2 + [z - 


(2n - l)h - s'cos(x 


- «>] 2 } 


+ Jz - (2n - l)h - s'cos(x - a)J(cos a) 


w'Kj 

S’ 

' co'^y 2 +| 

z - (2n - l)h - s'cos(x - 

■ otj 

i 5 } 

t 

Iy 2 + 

•[- 

(2n - l)h - s'cos(x - a)] 

2 


Wy 2 + [z - 


(2n - l)h - s'cos(x 


- a)] 2 > ^ 


(C12) 


A second exponential Fourier transformation on z - 2nh + s cos(x - a) with variable of 
transformation g and on z - (2n - l)h - s'cos(x - a) with variable of transformation 
g' gives by use of formula 43, page 56, and formula 43, page 112, of reference 15 the 
transforms corresponding to expressions (Cll) and (C12), 


dG n = (- l) n Aw 0 ds 


-i(sin Q!)a>e 


-lyl^+w 


^ -|y|yw 2 +g 2 

_ i(cos opge 


^g 2 + a> 2 


^co 2 + g 2 

(|y| >0; |w| >0) (C13) 


d^n = (-l) n AwQ ds' 


, -|y|Vg'“+w' 


' 2 +a >' 2 -I y I \/uj , 2 +g ’ 2 

i(sin oi)a)'e 1 J 1 1,0 * i(cos aQg’e 1 1 


\Jg ’ 2 + co 


,2 


\f7 2 + g' 2 


(|y| > 0; |co’| >0) (C 14) 
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The boundary condition (CIO) transforms to 

G±l §2 = 0 

and the Laplace equation — ^ = 0 to 

9x 2 9y 2 9 z2 


(Cl 5) 


9 2 G 


9y : 


2 = (^ 2 + g 2 )G 


(C16) 


A solution having the required symmetry in y is 

G(co,y,g) = B(o>,g)cosh(y||co 2 + g 2 ) (C17) 

and if dG n>s = B s (oj,g)cosh^y'Jc<j 2 + g 2 ^ is the influence of the slotted boundaries such that 
(dG n + dG n?s ) satisfies the boundary condition (eq. (Cl 5)), it follows that at the boundary 


v = - 
y 2 ’ 


B s(^,g) 


cosh^|)Ja> 2 + g 2 ) + l ||o> 2 + g 2 sinh^w 2 + g 2 j 


n . -|\|cu 2 +g 2 

= (-1) AwQds ie (o> sin a + g cos or) | 


Aj co 2 + g 2 


- l\ 


Similarly, if dGn,s * Bg cosh(y ^co’ 2 + g' 2 ), then 
B>^[cosh(|^7^) + Z\ju>' 2 + g' 2 sinh^^’2 + g’2^ 


O 1 * 

= (-1) Awq ds’ie (g'cos a - co’sin a)( 


|^’ 2 +g’ 2 


V 2 + g’ 2 


Satisfaction of the boundary condition at the boundary y = - | yields exactly the same 
equations. It therefore follows that 


B s (co,g) = 


/ xn -|\Jw 2 +g 2 

(-1) Aw Q ds ie * (co sin o?+ g cos 


”> f.r~r~Ty - 

w •>“ ± g" / 


cosh + l \jo> 2 + g 2 sinh ^u> 2 + g 2 ^ 


(|o)|>0) (CIS) 
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-|j)/w' 2 +g' 2 / j \ 

(-1) Aw 0 ds'ie 2 (g'cos a - co'sin a) . = = - £ 

, ■ )** + r - 2 ; (Kl>o) (C19) 

coshfl^co'S + g>2 J + i\J W '2 + g'2 sinhf^Va)' 2 + g'2j 

The transforms of the slotted-boundary interference potentials corresponding to the 
influence of the slotted boundaries on the differential disturbance potentials d(f> n and 
d(p' n are now 

n — D //.» « rr2^ 


dG n>s = B s (o>,g)coshlyVo^ + g 


2 4. or * 2 1 


J n,s = Bs(o)’,g')coshlyVw^ + g 


where B s and Bg are given by equations (C 18) and (C 19), respectively. The corre- 
sponding differential interference potentials are 


(-l) n Aw 0 ds 

d K s — llm A 1 + dw 

n > s 4 7 r 2 € —0 I I 11 


i|w[x-s sin(x-a)]+g[z-2nh+s cos(x-a^ ^ 2 ^ 2 +g 2 


(cu sin a + g cos a) / 


ta 2 + g2 


coshfe^) 2 + g 2 j + l^u) 2 + g 2 sinh(|^u; 2 + g 2 


cosh y vw 2 + g 2 ) dg 


( 6 >0) (C20) 


(-l) n Aw n ds' 

dtf)' = lim i 

n ’ s 4 tt 2 6-0 


i f [x- s ’ sin (x - «)] +g ’ [z - (2n- 1) h- s ' cos(x- a?)j} 


b^2+g^2 \ f I — o 9 ) 

e L (g’cos a - ui’sin at) — — - Mcoshlyya)’* + g’^J 

x Vj>'2 + g’2 j \ _ j dg , 

cosh^|^lu)'2 + g’2^ + j ^L'2 + g'2 sinh^^w' 2 + g ,2 j 


(e >0) 


(C21) 
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>-e r 

oo’ 

n-e r 

) + \ 

dco means 

\ dco + \ 

_oo 

■ 

-OO 


In equations (C20) and (C21), | \ + \ |dco means \ dco + \ dco. The limits 

evidently exist, at least for any finite value of l. 

The velocities in the x-direction on the vertical center plane y = 0 are now 


8 K,s) (-D^p ds 

8x 4,2 


-oo - oo 


i(^[ x-s sin (X" a )J + g[ z "2nh+ s cos(x-Q!)j^ - 


(co sin a + g cos a)(wl - = 


\ yu>2 + g2 j 


coshfWw^ + g2| + i y w 2 + g2 sinhfe^ + g 2 


dco dg 


(C22) 


9 ( d0 n,s) _ (-^Awp ds' 
9x 4 ^2 




V _ CO V-J 


i (u '[x- s ' sin(x~ a?)] +g ’[z - (2n- l)h- s ’ cos (x- Q!)]| 


x 


e-t^'Veos a . w . sin a) / u! . gl \ 

' — V dw , dg , 

cosh ^ \|co’2 + g'2j + £ ^'2 + gi2 s y* ^,>2 + 


(C23) 


The expressions for the velocities in the z-direction on the vertical plane are the 
same as equations (C22) and (C23) except that the quantities col — w and 

,^2Ti2 

co 7 — are reDlaced bv el. - — - and e-7 — H — . resDectivelv. To 

\Jco'2 + g '2 \/co'2 + g'2 

obtain the interference velocities due to the influence of the slotted side walls on a vortex 

cylinder cell, these expressions must be integrated over s and s', where s and s' 

run from 0 to - sec(x - a). The integrations of expressions (C22) and (C23) on s give 
z 
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,n) g 9(0 n + ^h) s _ (-D n Aw 0 

K 47T 2 i 


i[wx+g(z-2nh)J i jwjx- |tan(x-ar)J+g(z-2nh+|^ 


sin(x - of) - g cos(x - Qf) 


_ OO 


_k\| w 2 +g 2 

e ^ (w sin a + g cos a) /coZ - 


u 2 + g 2 


X 7 r - — — 7 v ■ - f d g 

coshfj^co 2 + g 2 j + l]ju) 2 + g 2 sinh^a) 2 + g 2 j 


i^ t x+g’|z-(2n<l)h^ i jupjx- |tan(x-a)|+g , |z-(2 n -l)h- ^ j 


a)’sin(x - Qf) + g f cos(x - o?) 


-vyu)’*+g'* . \ \ 

e 2 (g’cos a - aPsin a)(oo'l — w A ' 

\ V ^’ 2 + g ’7 . , , 

X — r - — - — ■ f daP dg’ 

coshf|\JuP 2 + g ,2 j + Z\|uP 2 + g ’ 2 sinhf]^aP 2 + g ^ 2 j / 


An obviously similar equation holds for the interference velocity in the z-direction. 

Equation (C24) and the similar equation for the velocity in the z-direction must be 
expanded into form suitable for computation as follows: 

r° P 00 j~[u)^sin(x - ar)sin a+ g 2 cos(x - a^cos aJ^sin(wx)cos|g(z - 2nh)j 


n|a> x - | tan(x - a)^cos g(z - 2nh + 0 j+ u)g£sin(X - nr)cos a 


-5Cv ; 


+ cos(x - ar)sin aj/cos(cox)sinJg(z - 2nhjj 


8 (»n+<frn) s (- l) n Aw Q V jj 2 + g 2 ] 

811 " 2 J 0 Jo 


Cjx - | tantx - <j)jjsin|g^ - 2nh + 

u> 2 sin 2 (x - a) - g 2 cos 2 (x - a) 


-[ui 2 sin(x - a)sin a+ g 2 cos(x - t*)cos ajj|sin(u>x)cos jgjz - (2n - l)hjj - sin|u>|x - | tan(x - a)J|cos|g|z - (2n - l)h - |Jj 
+ wg[sin(x - «)cos a+ cos(x - at)sin aj^cos(u)x)sin ^|z - (2n - l)hjj - cos {^[ x ■ I tan k ‘ " * 2n _ 1)h " ij}j 


ot^sin^ - a) - g2 cos 2( x _ 
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O^n+Qs (- l) n Awp 
bz 7T^ 


-|^> 2 +g2 

e z (g* - 




T?) 


cosh (l^ 2 + g2 ) + 1 + g2 sinh ^|^" + g 2 ^ 


[w 2 sin(x - a) sin a+ g 2 cos(x - a)cos arj|cos(wx)sin[g(z - 2nh)] 

- cos ju[x - | tan( x - a)jj sinjg(z - 2nh + + cugjjsin( x - a)cos a 

+■ cos (x - a) sin aj jsin(wx)cosjg(z - 2nh)j 

- sinjwjx - | tan(x - or)jjcosjg^z - 2nh 4 j 


te 2 sin 2 ( x - a) - g 2 cos 2 (x - a) 


-jaj2sin(x - a)sin a+ g 2 cos(x - or)cos a| |cos(wx)sin|gjz - (2n - l)hjj - cos jwjx - | tan( x - a)jjsin^Jz - (2n - l)h - |jjj 
+ a>g|sin(x - a)cos a 4 cos(x - a)sin aj^sin(u)x)cos|g|z - (2n - l)hjj - sin jw|x - | tan( x - or)j j cos {f[ z - (2n ' 1)h - 1]}) 


w 2 sin 2 (x - a) - g 2 cos 2 (x - «) 


d»><«g (C26) 


By equations (3) and (4) 

6 - u C 

R»x ~ w Q a 

and u 

. _ _v_C 

°R,z Wq A 


(C27) 

(C28) 


where u and v are interference velocities in the x- and z-directions, respectively, 
and C = hb is the cross-sectional area of the test section. Then with p = cuh and 
q = gh, the contributions (A6 R and (A6 RjZ ) ng to 6 R x and 6 R>Z , respectively, 

of the slotted -side-wall interference at y = 0, z = 0 on a vortex- cylinder cell are 
given by 


( A6 «.x)„ 


_(-l) K b 
> »2 h 


'>00 


p2sin(x - a)sin a + q 2 cos(x - <x)cos aj|sin^p ^cos(2nq) 


e -\X^Ll p \ 
( h & + 

- sin^pj^ - |tan(x - a)Jjcosjq^2n - |jj j - pq[sin( x - or)cos a 
4 cos(x - «)sin oj^os^p ^sin(2nq) - cos<jpJ| - | tan( x - a)jjsinjq^2n - |jj J 

0 w 

o cosh (i h^ 2 + + h^ p2 + q2 sinh (z ^ p2 + ^ 

p2sin 2 ( x - a) - q 2 cos 2 ( X - a) 


|j) 2 sin(x - ar)sin or 4 q 2 cos( X - <*)cos a)J ^sin^p ^Jcosjq(2n - l)j - sin cosjq ^2n - j 4 pqjsin( x - or)cos a 

4 cos(x - of)sin «]|cos^p ^sin[q(2n - lj] - cosjpj^ ‘ \ tan (x - oi)jjsinjq^2n - |jjj 


dp dq 


(C29) 
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{ A! 'R,z) na = 


Jo 




h &+q 2 / 


0 cosh(i |\|p 2 + q 2 ] + ^l(p 2 + q 2 sinh^ |Jp 2 + q^L 


Jp 2 sin(x - a)sin a+ q 2 cos( X - a)cos a J^ os |p|^ " ^ tan (x - «)j^ sinjq^2n - ^ 
- cos^p j^sin(2nq)^ + pq£sin( X - a)cos a + cos(x - a)sin aj^sin^p ^cos(2nq) 

' sin {*jj ' 1 tan(x ' «)jj cos [q( 2n - 1)] 


p 2 sin 2 (x - a) - q 2 cos 2 (x - a) 


[j) 2 sin(x - a)sin a + q 2 cos(x - a)cos ajjcos ( p h) sin [ q(2n " 1 )] ' ^ os {p[k \ tan h " 0( )jJ sin jq( 2n ‘ + pq|sin(x - a)cos « 

+ cos(x - «)sin aj|sin^p 2jcos[q{2n - 1)| - sin|pj| “ | tan (* - a)|jcosjq^2n - 


p 2 sin 2 ( x - a) - q 2 cos 2 ( x - a) 


dp dq 


(C30) 


For p and q both zero, the integrands in equations (C29) and (C30) are zero. For q 
approaching ±p tan( x - or), the limits exist, and since the functions are continuous, no 
difficulty should arise in integrating across the singularity provided proper attention is 
given to accuracy in the neighborhood of the singular point. 

For the condition p sin( x - a) approaching ±q cos( x - a), that is, q = ±tan( X - a), 
the quantities in braces in equations (C29) and (C30) become indeterminate and must be 
replaced as follows: 


For tan( x - «) positive and q/p approaching tan( x - a), the quantities in 
braces in equations (C29) and (C30), respectively, are replaced with 


q cos a ± p sin a 


4 cos( x - a) 


Cos^p ^ - 2nq^ - cosjj? 5 + (2n - l)qj 


+ sin p 


q cos a - p sin a 


2q cos( x - a) 


|cos£(2n - l)qj - cos(2nq)j 


(C31) 


and 


q cos a + p sin a 


4 cos( x - a) 


Qos^p 5 - 2nqj + cosjp j + (2n - l)qjj 


C0S ( P C 


q cos a - p sin a 


2q cos(x - a) 


|sin(2nq) - sin|^2n - l)qj 


(C32) 
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For tan(x - a) negative and q/p approaching -tan(x - o') the quantities in 
braces in equations (C29) and (C30), respectively, are replaced with 


~q cos a - p sin «] [ eo8 T 5 _ (2n . l)q] - cosfp £ + 2nq\ 
4 cos(x - a) l l “ J \ h ) 




0 


j 


(C33) 


and 


q cos a. - p sm a i 
4 cos(x - Qf) 



+ c°s(p| 


geosaH-pslnair (2 , . sin r (2n . 1)q j 

2q cos(x - a) L J 


(C34) 


For tan(x - a) - 0 and q = 0 the quantities in braces in equations (C29) and 
(C30) are to be taken as zero for any value of p. The total contributions of the slotted 
side walls to the interference factors are obtained by summing expressions (C29) and 
(C30) over all values of n: 


( 6 R,x) s - Z( A6 R,x)ns 

n 


(C35) 


( 6 R,z) s "Z( A6R ’ z )ns 

n 


(C36) 


To these contributions of the slotted side walls must now be added those of the 
images in the closed top and in the open lower boundary. For this purpose, expres- 
sions (C5) and (C7) must be integrated over s and s’, respectively. 


For the integration of expression (C5) over the original disturbance vortex sheet 
extending from the rotor disk to the free boundary (see fig. 15) let 


= 


Awq sin a 
2v cos(x - a) 


(C37) 


Then 


z 0 = s cos(x - a) 

s= fo 

cos(x - a) 


(C38) 

(C39) 
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s sin(x - a) = Zq tan(x - a) 

and the potential corresponding to the first term in equation (C5) is 

ph/2 

<t> i = 


[x - z 0 tan(x - a)]dz 0 


0 


j[x - z q tan (x - a)] 2 + y 2 + (z + z 0 y 
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(C40) 


(C41) 


where h is the test- section height. 

For the integration of the second term in expression (C5), let 

Awq cos a 


M 2 2t7 cos(x - a) 


(C42) 


The integral is then 

$2 = M 2 

With 


'h/2 


(z + z 0 )dz Q 


{jx - z Q tan(x - a)] 2 + y 2 + (z + z 0 )^ 


(C43) 


Zq = s' cos(x - a) 


(C44) 


similar expressions are obtained for the terms of the potential of the image in the open 
boundary. (See fig. 15.) By integration of expression (C7) 


= - Mj 


ph/2 


[x - Zq tan(x - opjdzfr 


|[x - Zq tan(x - a)] 2 + y 2 + (z + h - z^ 


3/2 


(C45) 


and 


sh/2 


<$>2 ~ m 2 


( z + h - z^dz^ 


o (f " z o tan ^ x " “i ] 2 + y 2 + ( z + h - z o)^j 


(C46) 


The velocity potential of the vortex cylinder cell nominally located at (0,0,0) is thus 


4>o ~ + $1 + $2 + ^2 (C47) 

The potential (-l) n $ n of a vortex cylinder cell with vorticity of the same sense nomi- 
nally located at (0,0, 2nh) is obtained from equation (C47) by replacing z in the expres- 
sions for <p 2 , and ^ with z - 2nh. The integrations to obtain cp 0 an d 

the corresponding streamwise and vertical velocities follow. 
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where 


h = 


>h/2 


^0 * I 1 + *2 + *3 + *4 


+ MgZ^dZQ 


"A 


jx 2 + y 2 + z 2 + 2[z - x tan(x - a)jz 0 + [j- + tan 2 (x - a)] 2 ^ ^ 


r> 


h - 


h/2 


[m 2 - Mi tan(x - a)]z 0 dz Q 


r ^ 3/2 

<x 2 + y2 + z 2 + 2 |z - x tan(x - a)jzg + Jl + tan 2 (x - a )] z o j 


> (C48) 


h = 


r h/2 1 

^M 2 z + M 2 h - Mjx 

) dz 0 

Jo i 

|x 2 + y 2 + (z + h) 2 - 2| 

[z + h + x tan(x - a)j 

z o + 

1 + tan 2 (x - cn)j 

4 2 


Oh/2 


14 = 


Jo 


[Mj tan(x - a) - M 2 ]z^ dz^ 


|x 2 + y 2 + (z + h) 2 - 2jz + h + x tan(x - a)jz^ + Jl + tan 2 (x - a )] z o^| 
Then for x, y, and z not all zero in Ij and I 2 , 

2 ^ 2(1 + tan 2 (x - «)] z 0 + 2jz - x tan(x - £y)^J ^M^x + M 2 z ^ 


3/2 


4(x 2 + y 2 + z 2 ) [l + tan 2 (x - a)] - 4 l[z - x tan(x - + y 2 + z 2 + 2^z - x tan(x - «)J z g + [l + tan 2 (x - o^ZQ^j 

^M x x + M 2 z^ 1 + tan 2 ( x “ «)] | + [ z ~ x tan (x “ a)]j 

|(x 2 + y 2 + z 2 )[l + tan 2 (x - a)J - [z - x tan(x - of)] j^|x 2 + y 2 + z 2 + 2[z - x tan(x - o)| ^ + Jj. + tan 2 (x - a)] 

(M jx + M 2 z^ j^z - x tan(x - a)] 

:2 + y2 + B 2)[! + tan 2 (x - g;)J - jjz - x tan(v - + y 2 + z 2 ^ 


h/2 
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l2 = - 


*1 

M 2 - Mj tan(x - arj] ^[z - 

x tan( x - a) 

z g + 2(x 2 + y 2 + z 2 )} 





j^(x 2 + y 2 + z 2 ) jj. + tan 2 (x - 

«)] 

- 4 

z - x tan(x - 

ajj^j jx 2 + y 2 + z 2 

+ 2 

z - x tan(x - <*) 

N 

O 

+ 

[l + tan 2 ( x - o;)j 

1 2 ^ 

z 0 J 

U2 


|h/2 


— JO 


[m 2 - M 1 tan(x - of j] ([z - x tan( x - c*)]| + (x 2 + y 2 + z 2 jj 


|(x 2 + y 2 + z 2 )[l + tan 2 (x - a)] - £z - x tan( x - Qf)J^|x 2 + y 2 + z 2 + 2[z - x tan( x - a)] ^ + |l + tan 2 ( X - ajj 

^2 - tan( x - a)J(x 2 + y 2 + z 2 ) 

|(x 2 + y 2 + z 2 )[l + tan 2 ( x - a)] - [z - x tan( x - or)] 2 | (x 2 + y 2 + z 2 ) 


1/2 


13 = 


2 + tan 2 ( x - o)|zq - 2^z + h + x tan( X - ajjj 


h/2 


|4 x 2 + y 2 + (z + h)^][l + tan 2 (x - a)] - 4|z + h + x tan(x - ajj^||x 2 + y 2 + (z + h) 2 - 2^z + h + x tan( x - aj] Z q + |l + tan 2 ( x - «)] z o^| 



( 

m 2 

z + M 2 h - M]X^ | 

[l + tan 2 (x - at)] 

| - [z + h + x tan(x - a)] j 

f 

||x 2 + y 2 + (z + h)^ 

* 1 

1 

1 

+ 

j j 

l-l 

|^z + h + x tan(x - 

■ a )]^jx 2 + y 2 ■ 

4- (z + h) 2 - 2 

|z + h + x tan(x + «)| 

IH 

[l + tan 2 (x - <*)] 

ih 2 l 

r 



^M 2 z + M 2 h - 

M i x )| 

[z + h + x tan(x - 

->] 


j^c 2 + y 2 + (z + h) 2 ] 

[l + tan 2 (x - 


|z + h + x tan(x - Of)j 

i?i 

[x 2 + y 2 + (z + h) 2 ] 

r 


I4 - ■ 


2[Mj tan(x - «) - Mj] 

-2^z + h + x tan(x - o)]zq + 2 [ x2 + Y 2 + ( z + h) 2 ] 

1 

j^x 2 + y 2 + (z + h) 2 ]^ + tan 2 (x - a) 

- 4|z + h + x tan(x - 

jx 2 + y2 + ( z + h)2 _ 2 jz + h + x tan(x - 

“)] z 0 + 

_1 + tan 2 (x - Of) 

.ji 

z 0j 

1/2 


n h/2 


|mi tan(x - a) - Mgjjjz + h + x tan(x - <a)j ^ - [x 2 + y2 + (z + h)^jj 


|^ x 2 + y2 + ( z + h) 2 ][l + tan 2 (x - a)J - [je + h + x tan(x - a )]^j|x 2 + y2 + (z + h) 2 - 2jz + h + x tan(x - a)] jj + [l + tan 2 (x ~ <*)] 
[mi tan(x - a) - M^jjx 2 + y 2 + (z + h) 2 ] 


|[x 2 + y 2 + (z + h) 2 ][l + tan 2 (x - a)j - [z + h + x tan(x - a |)^| (x 2 + y 2 + (z + h ) 2 ] 
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Then 


$0 * *1 + l 2 + *3 + l 4 



[l + tan 2 ( x - a)]| + [z - x tan( x - <*)] 

} - [m 2 - ] 

Mjtanfe-ajjl 

[z - x tan(x - a)]| 

+ (x 2 + y 2 + z 2 )j 

B 

jc 2 + y 2 ) + 2xz tan( x - a) + (y 2 + z^tan 2 ( x - cr)Jjjx - | tan( x - , 


y72 


[m 2 ^ M! tan(x - <*)) (x 2 + y 2 + z 2 ) - (M t x + M 2 z)[z - x tan( x - a ) ] 
[(x 2 + yty + 2xz tan(x + (y 2 + z2)tan 2 ( x - a) J(x 2 + y 2 + z 2y 2 


|m 2 (z + h) - MjxJ 

jjl + tan 2 (x - - £z + h + x tan(x - ajjj - |m 2 - tan( X - o;)j|jz + h + x tan( x - a)]| - [ 

x 2 + y 2 + (z + h) 2 jj 

§ 

c 2 + y 2 ) - 2x(z + h)tan(x - a) + [y 2 + (z + h) 2 ]tan 2 ( x - x - 1 tan( x - a) 


^1/2 


t [ M 1 tan (x - a) - M 2 ][x 2 + y 2 + (z + h) 2 ] - [Mjx - M 2 (z + hj][z + h + x tan( x - a)] 
{(x 2 + y 2 ) - 2x(z + h)tan(x - a) + [y 2 + (z + h)^jtan 2 ( x - a)j(x 2 + y 2 + (z + h) 2 ] 1/2 


and the corresponding vertical velocity is 


9z 


Mjjx + {2z + |)tan(x - a) - M 2 x - | tan( x - a)Jtan( x - ot) 


[(x 2 + y 2 ) + 2xz tan(x - a) + (y 2 + z^)tan 2 (x - a)] jjx - | tan (x - <*)J + y 2 + (z + 1) ^ 


1 ( z + |) x + | z tan (x - <*) + (y 2 + z 2 ) tan(x - a) + M 2 fr| - z)x tan(* - a) + | z tan 2 ( x - a) - x 2 - y‘ 


x ^2x tan(x - a) + 2z tan 2 (x - a)]jjx - | tan( x - Of)j + y 2 + (z + |^j + [x 2 + y 2 + 2xz tan(x - ot) + (y 2 + z 2 )tan 2 ( x - Q?)J( Z 

r- _ o ^i 2 / 2 


[(x 2 + y 2 ) + 2xz tan(x - a) + (y 2 + z 2 )tan 2 (x - - | tan (x - «)| + y 2 + (z + J* 


M 2 x tan( x - a) - x + 2z tan( X - a)J 


[(x 2 + y 2 ) + 2xz tan(x - ot) + (y 2 + z^)tan 2 (x - ot) j(x 2 + y 2 + z 2 )^ 2 

j[-xz - (y 2 + z 2 )tan(x - a)] + M 2 [x 2 + y 2 + xz tan( X - jj 2 * tan (x - ot) + 2z tan 2 ( x - Of)] [x 2 + y 2 + z 2 J 
[(x 2 + y 2 ) + 2xz tan(x - ot) + (y 2 + z 2 )tan 2 ( x - o?)j^ 


[(x 2 + y 2 ) + 2xz tan(x - a) + (y 2 + z^|tan 2 (x - a)j (x 2 + y 2 + z 2 ) 


73/2 


(Equation continued on next page) 
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x -(2z + | h) tan(x - a) + M 2 || tan(x - ot) - x tan( x - a) 
x 2 + y 2 - 2x(z + h)tan( x - a) + [y 2 + (z + h) 2 ]tan 2 ( x - Qf)J jjx - | tan( x - a)j + y 2 + 

L x + |)x - (z + h)(z + |)tan( x - a) - y 2 tan( X - a)] + M 2 f(z + h)£ tan 2 ( x - a) - (z + 3 |)x tan( x - a) + x 2 + y 2 


x f[-2x tan( x - a) + 2(z + h)tan 2 ( x - fl)]/ x - | tan( x - a) J + y 2 + + | 


+ |(x 2 + y 2 ) - 2x(z + h)tan( x - of) + [y 2 + (z + h)^]tan 2 (x - a^z + 


j(x 2 + y 2 ) - 2x(z + h)tan(x - a) + [y 2 + (z + h) 2 ]tan 2 ( x - a)j jjx - | tan( x - «) j + y 2 + (z + | 

M 2 x tan( x - oi) - Mjjx - 2(z + h)tan( x - a)j 

j(x 2 + y 2 ) - 2x(z + h)tan( x - a) + [y 2 + (z + h) 2 Jtan 2 ( x - of)l |x 2 + y 2 + (z + h) 2 J 


M^|[y 2 + (z + h) 2 ]tan( X - a) - (z + h)xj + M 2 |jz + h)x tan( x - a) - x 2 - y 2 jV^-2x tan( X - Of) + 2(z + h)tan 2 ( x - a)j 


(x 2 + y 2 + (z + h) 2 J + (z + h)^x 2 + y 2 ) - 2x(z + h)tan( x - at) + [ y 2 + (z + h) 2 ]tan 2 ( x - Qf jjj 

j(x 2 + y 2 ) - 2x(z + h)tan( x - a) + [y 2 + (z + h) 2 ]tan 2 ( x - a)j |x 2 + y 2 + (z + h)^] 


The induced velocity in the direction of the tunnel stream is 

"a^oj Mi( Z + |)-M 2 (z-|)tan( x - a ) 

[y 2 + (y 2 + Z 2 )tan 2 ( x - Of)] ^ tan 2 (x - Of) + y 2 + ^z + 

{m^| z + y 2 + z 2 jtan( X - a) + M 2 || z tan 2 (x - a) - y 2 


x ^2z — - tan 2 (x - a) + y 2 + ^z + tan( x - a) - |[y 2 + (y 2 + z 2 )tan 2 ( x - a)jtan( x - a) 
[y 2 + (y 2 + z 2 )tan 2 ( x - rt)] 2 j- tan 2 ( x - a) + y 2 + (z + 


(Equation continued on next page) 
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For the condition y/h - 0 with tan(x - a) * 0 the right-hand side of equation (C49) 
must be replaced with 

b 1 sin a f z/h 2 {h + ^) (h * *) 1 c 


h 2v cos(x - a) 2 


h 2, 


£+ 1 


Similarly, provided not all x, y, and z values are zero, 


(C50) 


[j + (2 ~ + Stan( X - a)] sin a - \ tan( x - a) tan( x - «)cos a 

j_hb^o/x y z\ b 1 L h v h 2J _ J =! — 

w 0 A 3z \h’ h’ hj h 2rr cos(x - a) f « /irt2 _ . . f /v ,2 /K vi| . . IE 1 _,1 2 . /l 


!k) 2 + {if + 2 K K tan( * ' a> + [(h) 2 + (k) 2 J-” 2( x - 4 {ft ' tan(x ' + (hf + (h + 5) 


(Equation continued on next page) 
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(Equation continued on next page) 
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(C51) 


Singularities in this equation with ^ = 0, %■= -2n (n * 0) can be avoided by taking 

n n 

- ^ 2n tan(v - a) and - # (1 - 2n)tan(x - a). This is easily done except for the condi- 
h h 

tion ^ = 0, tan(x - a) = 0. 

Reference to equations (C37) to (C47) shows that for the condition - = 0, 

h 

tan(x - a) = 0, 



and 



y=o 

tan(x-Q!)=0 



Aw 0 cos a r 1 1 h 
2,co S (x-<*)^ + h)2]3/2 
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so that 


_1_ hb 
w 0 A 


8 ^ 0 



tan(x-a)=0 


b cos a? 
h 277 cos(x - a 



1 


21 


3/2 




z/h 



J 


(C52) 


For the interference corresponding to the location of a vortex cell at (0,0,0), the poten- 
tial <t>Q has to be modified by subtraction of the potential of a vortex cylinder extending 
from s = 0 to s approaching «. Thus for the interference potential in the X, Z 
plane due to the original vortex cylinder cell, Ij and I2 of equation (C48) must be 
replaced with 


M 


^M]X + M2z)[l + tan 2 (x - a)] 


y=o 

n =0 


j(x 2 + z 2 )|~l + tan 2 (x * Qf)| - [z - x tan(x - + tafl2 (x - <*)] ^ 

(M x x + M 2 z) jjl + tan 2 (y - a)]| + [z - x tan( x - <*)]j 
|(x 2 + z 2 ) jjL + tan 2 (x - on)] - - x tan(x - + z 2 + 2jz - x tan( x - a)J| + [l + tan 2 ( x - 


Ny=o = ~ 


y=o 

n =0 


Jm 2 - Mj tan(x - q)j jjz - x tan(x - qjj j + x 2 + z 2 


2 + z2 )t + tan 2 ( x - q)J - [z - x tan(x - ajj^||x 2 + z 2 + 2 jz - x tan(x - ajj| + [l + tan 2 (x - af| ^ 
[M2 - tan(x - a)jj~z - x tan(x - ar)j 
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j(x 2 + z 2 ) jl + tan 2 (x - q)j - |z - x tan(x - «jj j 1 [l + tan 2 (x - a)] 


1/2 


These quantities are obtained from Ij and I2 by changing the limits of integration to 
h/2 and It follows from these equations and the expressions for I3 and I4 that 
the velocity in the z-direction on the X-axis due to the semi-infinite vortex cylinder 
extending below the position of the open boundary being subtracted from the image in the 
open boundary is 


61 



This equation contains singularities for x = 0. Since these singularities result from the 
integration of the part of the vortex cylinder arising from the helicopter rotor, it is con- 
venient to start again from equation (C5) and write (with negative sign because of sub- 
traction of this part) with y = 0 the potential at y = 0 of the part of the rotor wake 
outside the position of the open lower boundary as 
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= 3 



[m 2 - M x tan( x - a)]dz 0 
[l + tan 2 ( x - a)J 5 /, 2 z 0 3 



Consequently, performing the integrations and adding in the terms contributed by I 3 and 
I 4 at x = y = z = 0 yields 
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Similarly, 


r a , n 


^00 

i 8< y = 

'j 


mgmm 

■1 

fz 

h/2 t- 


dz n 


^372 + 3M l 


1 h/2 


[- y,,<, - . jV, 

[z 0 2 tan 2 ( x - a) + z Q 2 ] 


h/2 


z 0 2 tan(x - «)dz Q 
[z 0 2 tan 2 (x - a) + Zg 2 ] 
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n=0 


= 3 


h/2 


[Ml tan(x - a) - M 2 ]tan( x - a)dz Q 
[l + tan 2 (y - «)] 5/2 z 0 3 


dz n 


[l + tan 2 (x - «)] 3/2 z 0 3 


h/2 


and 


8 \wl 2 frl 

tan(x - ot) - M2]tan(x - a-) 

M X 

. Jx=0 h 2 I 

y=0 L 

[l + tan 2 (x - a)] 5 ^ 2 

[l + tan 2 (x - a)\ 3 J 


M^-SM^tan^-a) 


n=0 


^tan 2 (x-a) + ^ 


V2 


Jmj y tan(x - a) - M 2 y tan 2 ( x - a)j|-2h y tan 2 { x - a) + yj - y tan 2 ( x - aytan( x - a) 

r i 3//2 

[h 4 tan 4 ( x - «)j|y tan 2 (x ‘ “) + yj 

M 2 h tan( x - a) - Mjh [Mjh 2 tan( x - <x)][-2h tan( x - a)] 
h 2 tan 2 (x - a) 


h 5 tan^(x - a) 


In form for computation the velocities corresponding to the image in the open lower 
boundary with the subtraction of the part of the vortex cylinder outside the test section for 
x = y = z = n = 0 are 


1 hb 

p>vl 


_ b 1 

w o A 

BZ _ 

x~0 

h 2 it cos(x 


3 sin a - tan(y - oQcos a 


y=0 

z=0 

n=0 


fsin a - tan(y - ar)cos J[2 tan 2 (x - ot) + 4_ _ sin a 
[l + tan 2 (x - a)J 3/,2 tan(x - a) 


(C53) 
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1 hb| o<p 0 T w| b 1 J6|tan(x - a)sin a - cos o)tan(y - a) o ^ 0 ^ / \ 

= r« 7 ?(— -=j A 2 sin a . sin a - 3 tan(y - a)cos a 

0 A L ®xJx=0 h 2ir eos(x -an r , -i5/2 f ^ ' 'n3/2 r I TiT? 

y=0 L L 1 + tan (x ■ «)J L 1 + ta “ 2 (x - “)] ' [1 + tan 2 (x - a)] 1 ' 2 tan 2 ( x - ct) 

z-0 J - J ' 

n=0 

2[sin a - tan(x - a)cos a]|~l + 2 tan 2 ( x - a)] ^ tan( x - a )cos a - sin a 2 sin a 1 , . 

(l + tan 2 ( x - a)J 3,/2 tan 2 ( x - a) tan 2 ( x - a) tan 2 ( x - a) j 

For tan(x - a) = 0, the quantity in braces in equation (C53) must be replaced with 

cos a (C55) 

and the quantity in braces in equation (C54) must be replaced with 

- | sin a (C56) 

0 ( p 

In form for computation, the velocity f or which equation (C53) is the special- 

a Z Jy=0 


(C54) 


(C55) 


(C56) 


ization at x = 0 is 
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Singularities can be prevented by taking ^ * tan(y - a). The equations are now at hand 
for calculating the interference factors defined by equations (3) and (4) in certain parts 
of the field. In particular, the upwash interference factor 6 R z can be calculated at 
y = 0 and z = 0 and the streamwise interference factor 6p^ x at x = 0, y = 0, and 
z = 0. 


is 


For x * 0 the upwash interference factor on the tunnel center line y = 0, z = 0 


V*]y =0 

z=0 



(C58) 


where the quantities to be summed are calculated according to equations (C30) with 
(C32) and (C34), (C57), and (C51), and in equation (C51) ^ = 0 and ^ = -2n. 

For x = 0, y = 0, and z = 0, the upwash interference factor is 



(C59) 


where the quantities to be summed are calculated from equations (C30) with (C32) and 

(C34), (C53) with (C55), and (C51) with (C52), with ^ = 0, ^ = 0; and in equations (C51) 

and (C 52), ^ «= -2n. 
h 

The streamwise interference factor at x = 0, y = 0, and z = 0, obtained as for 
equation (C59) by use of equations (C29) with (C31) and (C33), (C49) with (C50), and (C54) 
with (C56), is 


6 "1 1 hb 

R > x Jx=0 w 0 A 


y=0 

z=0 



I <-»“ 


dct> 0 ( 0,0, -2n) 


n*0 


9x 


I[(*y*lL 0 < C6 °> 


Two further specialized applications of the theory of this appendix are required. 

The first is the upwash interference velocity on a lifting element centered between two 
infinite slotted side walls. This interference velocity can be obtained from equation (C22) 

CO £ . 

by replacing the factor coi — ■ with g l — = For convenience, n is 

\|u> 2 + g2 'Ico 2 + g 2 
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assumed to be zero. Also, to correspond with the wing theory, a is assumed to be zero 
and y is assumed to be n/2. Then, the upwash interference velocity is 



The range of integration on s is 0 to 00 and since for x * 0 


= lim 
0 K — 00 


e io>x _ e iw(x-K) e ia)x 


lim 


-icoK 
e 

ico 


where K is any large number, it follows that in consideration of the oddness in co and 
the evenness in g, 



^00 



sin cox 

+ lim 

sin coK 

10 

0 L 

K — 00 

CO 


cos(gz)e 




Substituting p for bio and q for bg and noting that the integrand involving 

lim ? . in coK contributes to the integral only for co approaching 0 allows this equation 
K — 00 " 

to be written as follows: 
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poo poo 


bv s_ 1 

Aw 0 " „2 


Sin ( P b) 


/ z\ -^P 2 -«J 2 
cosfqgje z 


qq 


i <1 


b Vp 2 + q 2 , 


cosh^p 2 + q 2 j+ ^ 'jp 2 + q 2 sinh^|\p 2 + q 2 j 


■ dp dq 


^0 

pOO O 0Q 


11m 

cosh ^ ~ q sinh ^ K-«= 


sin(pK) 


dp 


'-’0 


3in ( p § 


n 2 J _ n 2 


, * 
cos/q 




L i . \ \ 


l b vtz) 

, dpdq+ If 

! sinh^||p 2 + q 2 ^ 


c »( q g). 2 <i(q£-i) 

cosh 5 + - q sinh - 
2 b 2 


dq 


(C61) 


In this equation for p approaching 


!!±! 


approaches x/b. With p = 0 and q 


approaching 0, 


V? 


approaches 1. 


+ Q" 


The second specialized application of this theory involves simply a modification to 
permit calculation of the upwash interference factor 5j^ z for skew angle x 
approaching n/2 and with angle of attack a equal to zero. Because of some diffi- 
culty with nonuniformity of convergence for x approaching ir/ 2, it is necessary to go 
back to the elementary potentials. With a = 0, y = 0, and x = equations (C8) and 
(C9) become 


d^ n = (-1) 


n^O 

2ir 


z - 2nh 


and 


r. Aw n 


life - x) 2 + (z - 2nh) 2 ] 

z - (2n - l)h 

j\|(s’ - x)^ + jz - (2n - l)hj^ 


\ds’ 


To obtain (f> n and 0^, these expressions must be integrated on s and s' from 0 
to The results are 


and 


<Pn = (-1) 


_ Aw n 
4>' = (-i) n — 2 
V n v ; 2n 


n Aw 0 1 


2tt z - 2nh 
1 


1 + 


\fx2 + (z - 2nh)^ 


[z-(2n-l)h]<| p ~ ~ _ (2n _ 1)h -ji 
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The corresponding velocities in the z-direction at z = 0 are 

d<p n 

Tr=<-'> 


Aw q 

-1 

27T 

4n 2 h 2 


h 2 1 1 + \p v) ( x 2 + 4 „2 h 3) 3/2 


(C62) 


d<b' Aw nr < 

^n _ / i\n 01 -1 


8 Z 


-Mr 


2 . <( 2n - l) 2 h 2 ^ ^2 + (2n . 1 } 2 h 2 j [ x2 + (2n _ 1 } 2 h 2 ] 


1 + 


372' 


(C63) 


The upwash interference velocity corresponding to the images in the top and bottom walls 
is given by the sum of expression (C63) for n = 0 and expressions (C62) and (C63) 
summed over all other integer values of n. 

To obtain the upwash interference velocity due to the slotted side walls at 
y = z = a = 0 and y = ^, expressions like equations (C22) and (C23) except for the 

CO CO ^ 

replacement of 0 ol - ~—=z and w'Z + — with gZ - — — — and 

\C27P Vl0'2 + g'2 C27? 

g'Z — — — — , respectively, must be integrated on s and s' from 0 to 00 and 

\lcL» *2 + g'2 

summed over all 0, positive, and negative integer values of n. The integration of 

e -io>s gi ves 


r 


r iws ds = J-- lim e 


•1CUS 


1W s-°o la > 


so that 

8 (Vs + K,s) 


dZ 


(-l) n Aw Q 


00 


y=0 4ir 2 
z=0 

Q *=0 

7T 

X= 2 


^i(u>x- 2 nhg) + e i[^x-( 2 n-l)hg]| e - 2 ^ 2+ ^ 


kJ - 00 vJ - 0O 


g(X - lim C— )(g l 

\ 1W s _oo iw A 


V ^ 2 + g 7 


cosh^cu 2 + g^ + Z io 2 + g 2 sinh^cu 2 + g‘ 


do> dg 


From considerations of symmetry with substitution of p for u>h and q for gh, this 
equation can be written 
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The integrands involving lim of 

£_oo 

h 

only for p — 0. Moreover, lim 

S_oo 

h 



sinm pjj contribute to the integrals 



and it is easily shown that 




dp = 0. Therefore, 
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Since 


v where v is the total upwash interference velocity, taking the 


required sums of expressions (C62) , (C63) , and (C64) gives 





( C65) 
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APPLICATION OF CONFORMAL TRANSFORMATION IN THE TREFFTZ PLANE 
FOR CALCULATION OF LIFT INTERFERENCE IN SLOTTED TEST 
SECTION SYMMETRICAL ABOUT THE VERTICAL AXIS 

The disturbance due to the lifting wing is represented by vortices trailing horizon- 
tally downstream from the wing tips as in reference 13. The interference velocity is 
computed in the far- downstream (Trefftz) plane, where the theoretical flow is two- 
dimensional, by means of the conformal transformation of the rectangular boundary to a 
circle. This application is limited to slotted configurations symmetrical about the verti- 
cal Y-axis and to slot arrangements for which both X- and Y-axes cut only closed por- 
tions (panels) of the boundary. The number of slots is then even. The wing is assumed 
to be mounted with its span horizontal and centered on the Y-axis, but it may be located 
off the X-axis. 

Symbols peculiar to this appendix are as follows: 


aj,bj real constants to be determined 

F,Fg j,Fg 2 specifically defined functions of 6 , 
h height of test section 

imaginary part of 

i-fi 


0g j, and respectively 


K 

K’ 

k,k' 


m,q 


& 


quarter-period of Jacobian elliptic functions 
coperiod corresponding to K 
arguments of K 
number of panels or slots 

factor in complex potential determined by slot arrangement 
real part of 
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V 

V 

* w 

X,Y,Z 

A 

x,y 

z 

r 

&VV 

e 

? 

$ 

V 

e 

a 

X 

0 


tunnel stream velocity 

upwash interference velocity, positive upward 

complex potential 

coordinate axes 

distance along axes 

complex position in z-plane, x + iy 

circulation 

upwash interference factor for wings 
small positive quantity near zero 
complex position in £-plane, £ + it) 
real part of £ 
imaginary part of £ 
polar angle 

complex argument of exponential as defined in equation (D13)) 
constant of conformal transformation (see eq. (D2)) 
amplitudes corresponding to arguments of elliptic functions 


Subscripts: 


g number designating a slot edge 

n,j index numbers 

t designates position of vortex in right half-plane 

Horizontal bar over complex quantity indicates the conjugate. Inclusion between 
two vertical bars indicates absolute value. 
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For convenience the width of the rectangular test section is taken as unity. Its 
height is designated h. The conformal transformation which carries the rectangle in 
the z-plane into the unit circle in the £-plane is discussed in reference 13. For conven- 
ience and because of the extension required to the development, the procedure of this 
reference is somewhat modified. The transformation is given by 


sn dn 


? = ■ 




where with the width of the test section taken as unity, 

A = 2K = ?£ 


(Dl) 

1 


(D2) 


Since for every value of the (quarter) period K of the elliptic functions sn, dn, and cn 
there exists the corresponding coperiod K', the values of K and K' may be found 
from the relation 

= h (D3) 

K 

with the aid of tables of the complete elliptic integrals K(k 2 ) and the relations 

K’ =K(k' 2 ) (D4) 

with 

k’ 2 = 1 - k 2 (D5) 

Note that on the circle 

? = e ie (D6) 


where 9 is the polar angle. Therefore by equation (107) of reference 13 the slot edges 
6g on the circle correspond to the slot edges Zg on the rectangle through 



1 - cn K) 

1 + cn^VZgj 


(D7) 


From the relations of reference 17 (p. 503), it can be shown that with x real 

cn(x - iK’) = ik”* — . Therefore, on the lower boundary where 
sn x 

XZg = Xxg - = Xxg - iK', it follows that 



Lower boundary 



1 

| k 2 sn^XXg) - dn^AXg 

'j - 2ik sn( 

Xxg^dn^AXg^) 

4 . 

sn| 

f XX gy 

sn (Xx g ; 

1 k 2 sn 2 ^Xxg^ 

) 

+ dn 2 ^Xxg 

) 
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so that with positive Xg 


%] 


Lower boundary 


= " 2 Sln 


2k sn^Xxgj 

dn (Xx g ) 

k 2 sn 2 

(Xx g ) 

+ dn 2 

K) 


-Ism’ 1 


2k\|l - k 2 sin 2 <£ sin <pj 
(D8) 


for values of (p corresponding to values of Xxg in tables of elliptic functions with the 
given values of X and K(k^ . 

Similarly, at the right-hand side boundary Xz g = X^| + iyg'j = K + iXy g so that 



Right side 


1 - cn 
1 + cn 


( K + ix yg) 

(K + iXy g j 


1 + ik’ 


( x yg> k ’) 

dn(Xy g ,k') 


sn| 

1 _ it 


X 1JA 

dn 

(Xy g ,k') 


and since 

dn^Xy g ,k'^ + k’ 2 sn 2 ^Xyg,k’^ = 1 

it follows that 

e 4ight side ■ 2 sl "‘ - \ i^J sin ' 1 i 2k ' 'J 1 - k ’ 2sln2 f' sin V 

(D9) 


for values of <p corresponding to values of Xy g in tables of elliptic functions with 
the given values of X and K^k’ 2 ). Because of the symmetry about the Y-axis of both 
the rectangular-test- section configuration and the conformal transformation, equa- 
tions (D8) and (D9) make possible the location of all the slot edges in the transformed 
plane. Still to be determined are the positions of the vortices trailing from the wing tips. 

Let the center of the two-dimensional vortex in the right-hand part of the z-plane 
be located at ^xt,y^ so that the complex position is 

z t = x t + 


By equation (Dl) the location of this vortex in the £-plane is 


?t = k + ir ?t = sn 




(DIO) 


By reference to formulas for expansion of the elliptic functions of complex arguments, it 
may be shown that if z = x + iy, then 


75 


APPENDIX D 


and 


sn(x)dn(y) + i cn(x)dn(x)sn(y)cn(y) 
sn(z) = ~ ^ 

cn^(y) + k z sn^(x)sn^(y) 


(Dll) 


dn(z) _ [dn(x)cn(y)dn(y) - ik 2 sn(x)cn(x)sn(y)] [cn(x)cn(y) + i sn(x)dn(x)sn(y)dn(y)] 
Cn ^ [cn(x)cn(y)] 2 + j^sn(x)dn(x)sn(y)dn(y)] 2 


(D12) 


In equations (Dll) and (D12) the elliptic functions with argument x are taken from 

tables corresponding to the value of K (or k2) and those with argument y are taken 

from tables corresponding to the value of K' (or k'2). With ^ x* in the place of x 
X ' 2 L 

and — y^ in place of y, equations (Dll) and (D12) can be used with some complex 

algebra to calculate ^ and from equation (DIO) and give the position of the vortex 

in the right-hand ^-half-plane. The other vortex is then located at the point (-If^t)' If 

the span lies along the X-axis so that y^ is zero, equation (DIO) is easily used directly 

to calculate ^ or, even more simply in this case, equation (107) of reference 13 


gives = 


1 - cn 


N 


1 + cn(\xjj 


1/2 


t 


With the configuration in the £ -plane completely determined, the flow due to the two 
vortices in the presence of the slotted circular boundary can be investigated. With a 
defined by the equation 


C = 




(D13) 


where the complex a approaches the real variable 6 at the circular boundary, it is 

seen that * has singularities, which according to the theory of refer- 

'Jsin 6 - sin 0 g 

ence 13 correspond to slot edges, symmetrically arranged with respect to the Y-axis. 
Since this symmetry has been assumed and since the axes cut only panels, it follows that 
the basic equation of reference 13 applicable to this analysis is equation (15) of that ref- 
erence with q equal to the number of panels (or slots) m and p equal to i. How- 
ever, this equation is specialized to apply only for the span lying along the X-axis. The 
corresponding equation applicable for the span parallel to but not necessarily lying along 
the X-axis is 


dw 

d? 


m 


+1 


2 

p? y ja n cos(no) + b n sin(no)j 
n=0 L 


(? - ?t)(? + CtX 1 - « t )(l + « t ) ^(sin a - sin 0 g ) 


a/2 


(D14) 
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where £ = e ia , l is the conjugate of $, £ t is the position of the vortex in the right- 
hand half-plane, the a n and b n values are constants, w is the complex potential in 
the £-plane, and for this investigation p = i. At the circular boundary, where £ = e i0 , 


m , 

T +1 


dw 

dK 


ie" i0 ^ ja n cos(n0) + b n sin(nd)j 
n=0 


a=e (e ie - ? t )(e i0 + ?t)( e_i0 - ?t)( e ' i0 + Ct) 6 - sin 0 g ) 


1/2 



ie" i0 |a n cos(ne) + b n sin(n0)j 
n=0 


- C t | 2 | e _ie + ff^sin 9 - sin 0 g ) 


2 m 


A/2 


Then 


m.+i 

2 


le 


dw 

dz 


-i(77-0) ^ |a n cos(mr)cos(n6) - cos(nir)sin(n0)j 


n=0 


CfeTT- 6 


e i(ir- 6) _ ^ 2 e -i(ir-0) + ^| 2 n^sin 6 - sin 6 g) 1 / 2 


m+i 

2 


-ie^ ^ ja^ cos(n7r)cos(ne) - b n cos(n7r)sin(n0)j 


n=0 


.-ie 


ie 


-e -e + C 


2 na 1/2 

II ^sin 9 - sin 6^j 


and it is true as required by the symmetry that 


dw 

d£ 


dw 


provided 


l 2j = 

2j+l = 


(j = 1. 2, 3, 
(j = 1> 2, 3, 


. •) 
• •) 
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Consequently, for symmetry about the Y-axis with the axes cutting panels, equation (D14) 
becomes 


dw 

d? 


[09 

!(t ■ ► 1 


£ a 2n+l 

cos(2n+l)a+ ^ 

b 2n sin(2no) 

n=0 

n=l 



(? “ ?t)(? + ?t)(! - ?C t )(l + ^t)^( sin a “ sin %) 1/2 


(D15) 

1 


where m is the (even) number of slots and 


to f 


is the largest integer less than or equal 
1/m 


and [Kt + *)] is the largest integer less than or equal to + lj. At the cir- 
cular boundary, equation (D15) becomes 


le 


■ie 


dw 

d£ 


M 

1 a 2n+l cos(2n + 1)0 + ^ b 2 n sin(2n0) 


n=l 


0*0 


e i0 - ? t | 2 je“ i0 + ? t IT^sin 0 - sin 0 g J 


(D16) 


\l/2 


Since the real part of w is the potential and its imaginary part is the stream func- 
tion and since the panel is necessarily a streamline (stream function constant), equa- 
tion (19) of reference 13 is equivalent to 


f = 

^Panel ds 


0 


(D17) 


At the wall, £ = e iy and d£ = ie iy d0. Moreover, with ^ + ^t> ^ * s not difficult 

to show that 

|e 10 - ? t | 2 |e" 16 + $ t | 2 = (l + r t 2 ) 2 - 4| t 2 - 4?7 t (l + r t 2 )sin 0 + 4r t 2 sin 2 0 


where r^ = | ?t| = ^ 2 + With these substitutions the use of equation (D16) in equa- 
tion (D17) gives, since m is even, 


^ Panel 



rml 

M 

^ a 2n+l 
_n=0 

cos(2n 

+ 1) 0 + ^ b 2 n sin(2n0) 

n=l . 

d0 

”(l + r t 2 ) 2 - 4?, 2 - 4., t l 


o 9 

Isin 0 + 4r^ sin^0 

m 

11 

g=l\ 

jsin 0 - sin 0 g | 

if 2 


= 0 (D18) 
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The limits on the integral of equation (D18) are the values of 6g at the edges of any 
panel wholly in the right half-plane. Because of the symmetry about the Y-axis and of 
the fact that the Y-axis cuts panels, there are ^ - 1 independent equations (D18) linear 
in the constants aj and bj. The other two equations needed for determining the -j- + 1 
constants are obtained from equation (20) of reference 13, where the contour integral is 
taken about the vortex at ^ in the C- plane. 


Since 


and similarly 


e i(T - e“ ia 
sin a = 

2i 


i - r 1 

2i 


cos a = 


g + r 1 

2 


it is evident from equation (D15) that the residue of dw/d£ at ^ is 





a 2n+l 


2n+l 


+ ?,- <2n+1) " 

~2 



(ct+ CtX 1 + 


nkX 

g=l\ 2i 



Since the value of the contour integral is the product of 2ni by this residue, equa- 
tion (20) of reference 13 gives the two relations 



(D19) 
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where and ^ indicate real and imaginary parts, respectively, and the circulation r 
has been taken to be unity. 

The sign in equation (D19) is ambiguous because of the two branches of the half 
power of the n-product of complex quantities in the denominator. In practice either sign 
and either branch may be chosen and its correctness easily determined from physical 
considerations. Because the constant on the right-hand side of equation (D19) is the only 
nonzero constant appearing on the right of any of the equations in the simultaneous linear 
set determining the constants aj and bj, the choice of sign in this equation determines 
the signs of any velocities calculated. If these signs are incorrect, they need only be 
reversed. 

Equations (D18), (D19), and (D20) provide the ^ + 1 linear simultaneous equa- 
tions for determining the y+1 constants aj and bj. The calculations are best per- 
formed by a machine computing system, because for more than about four slots they 
become too complicated, voluminous, and sensitive to errors of various kinds to be suit- 
able for hand calculation. If the computing system permits use of complex variables, 
equations (D19) and (D20) may be used directly without the necessity of reduction to real 
functions of real quantities. It is only necessary to insert the limits on the summations 
and on the n-products along with 

*t = k + 
h = k ~ i71 t 

|S t | 2 = r t 2 =«t 2 + V 

and the values of Og. The coefficients of the aj and of the bj are then easily obtained. 
Thus, the coefficient of a 3 in equation (D19) is 
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The denominator must remain the same in all the expressions for the coefficients of the 
aj and bj in equations (D19) and (D20); that is, the same branch of the n-product must 
be used in calculating all the coefficients. 

The integrands of the integrals determining the coefficients of the aj and bj of 
equation (D18) have singularities at the limits of the nature of l/^e for e > 0 as e 
* approaches 0. If the numerical integration routine does not depend on the values of the 
integrand at the limits, correct values of the coefficients may possibly be obtained by 
simply ignoring the singularity, since the contribution to the integral of the region within 
e of either limit is of the order of \Je\ Otherwise, integrals of the form 


'%,2 / sin\ 

or j0 d0 


Vcosy 


g,l 


ne) 


e g,2 > e g,l S 

J = 1> 2, . . 


may be approximated by 


%,2 / sin\ 

I or j0 d0 
\cos/ 


0 g,l 
where 


no) 


'%,2-e / sin \ 

I or ] j 0 d0 
\ cos/ 


6 g,l +e 


no) 


+ 2 Ve 


( or W,2 

\cos/ 


( 

/sin\ 

or 

Vcos / 

j%,l 

F g,l 

(%,l)\|cos %, 1 


F(e) = 


+ r^ _ 4^2 _ 4rj t ^l + r t 2 ^sin 0 + -Ir^sin^ 


IT ^|sin 0 - sin Q 0 


? g 


F g, 2(%,2) " 


(l + r t 2) 2 - 4 ? t 2 - 4r? t (l + r^sin 0 g>2 + 4r t 2 sin 2 0 g>2 sin 0 gj2 - sin 0 g |) 

g^g2 


1/2 


,l( e g,l) = { 1 + r t 2 ) 2 - 4 ^t 2 - 4?? t( 1 + r t 2 ) sin e g,l + 4r t 2sin2 %,l ^ (|sin e gfl - sin 0 g | j 


J g= 

g^gl 


a /2 


and t (radians) is a small positive number. Depending on the width of the panel and 
the number of points in the integration routine, e may be taken to be of the order of 
0.0001 to 0.01. 
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Once the constants aj and bj have been obtained, the complex velocity dw/d£ 
can be calculated at any point £ in the £-plane by use of equation (D15), where 
cos(2n + l)a, sin 2na, and sin o are replaced by their expressions in terms of 
£ (= e i0 ). The complex velocity at the corresponding point in the z -plane is then 

For this investigation, the point of greatest interest in the z -plane is z = 0 and 
the corresponding point in the £-plane is £ = 0. At this point equation (D 15) simplifies # 
by letting £ approach 0 to 


dw 

d? 


21-1 

±i2 ^ /a m + b„ 

f +1 T +1 ' 


K=o 


(D21) 


where one of the constants is zero, depending on the value of m. The ambiguity in sign, 
which corresponds to the ambiguity in sign of the constants, can be settled by noting that 
for a lifting wing the velocity produced by the vortices at the center of the span is down- 
ward (that is, negative) at this point. It is also negative, and purely imaginary, at other 
points on the Y-axis, where the Y-axis bisects and is normal to the span (displaced to the 


dw 


Trefftz plane), so that for a lifting wing^— 


<=0 


is positive. 


X dw 


C=0 

Is. 


At z = £ = 0, ^ thus, the velocity in the z -plane at this point is -D ^ 

dz 2 " * a s_ 

This velocity is the total velocity due to the vortices in the presence of the slotted wal 
The interference velocity at z = 0 in the plane of the wing is half the difference between 
this velocity and that at z = 0 in the Trefftz plane due to the vortices in the infinite free 
field. Thus, the interference velocity is 

x t 


1 - 1 

- A ( 

) dw 

Jz=0 " 2 

2 d 

d? 


r=i 


*(*t : 


2 + y t 2' 


where the circulation r is taken positive for a lifting wing. For r * 1, it is only neces- 
sary to multiply the right-hand side of this equation by T, the negative sign being included 
if the wing is producing negative lift. Then 



r x t x 

fldw 

1 

Jz=0 " 2 

*( x t 2 + yt 2 ) 2 

IrdS 

S' s 

II 

l_2_ 


(D22) 


and since 


C L S = 


4r * t 

V 


it follows from equation (1) that 
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(D23) 


where for this appendix b = 1. Because it is total induced velocities that are computed 
and the interference velocity v is the difference between two such total induced veloc- 
* ities the computations must be carried out with second-order accuracy, that is, with such 
accuracy that the differences are significant. 


( 


83 


REFERENCES 


1. Wright, Ray H.; Ritchie, Virgil S.; and Pearson, Albin O.: Characteristics of the 

Langley 8-Foot Transonic Tunnel With Slotted Test Section. NACA Rep. 1389, 

1958. (Supersedes NACA RM L51H10 by Wright and Ritchie and RM L51K14 by 
Ritchie and Pearson.) > 

2. Heyson, Harry H.: Jet-Boundary Corrections for Lifting Rotors Centered in Rec- 

tangular Wind Tunnels. NASA TR R-71, 1960. t 

3. Heyson, Harry H.: Linearized Theory of Wind-Tunnel Jet-Boundary Corrections and 

Ground Effect for VTOL-STOL Aircraft. NASA TR R-124, 1962. 

4. Heyson, Harry H.; and Grunwald, Kalman J.: Wind-Tunnel Boundary Interference for 

V/STOL Testing. Conference on V/STOL and STOL Aircraft, NASA SP-116, 1966, 
pp. 409-434. 

5. Cook, Woodrow L.; and Hickey, David H.: Comparison of Wind-Tunnel and Flight- 

Test Aerodynamic Data in the Transition- Flight Speed Range for Five V/STOL 
Aircraft. Conference on V/STOL and STOL Aircraft, NASA SP-116, 1966, 
pp. 447-467. 

6. Anscombe, A.; and Williams, J.: Some Comments on High-Lift Testing in Wind Tun- 

nels With Particular Reference to Jet-Blowing Models. Rep. 63, AGARD, North 
Atlantic Treaty Organization (Paris), Aug. 1956. 

7. Rae, William H., Jr.: An Experimental Investigation of the Effect of Test Section 

Geometry on the Maximum Size Rotor That Can Be Tested in a Closed Throat Wind 
Tunnel. Paper No. 66-736, Am. Inst. Aeron. Astronaut., Sept. 1966. 

8. Theodorsen, Theodore: The Theory of Wind-Tunnel Wall Interference. NACA 

Rep. 410, 1931. 

9. Davis, Don D., Jr.; and Moore, Dewey: Analytical Study of Blockage- and Lift- 

Interference Corrections for Slotted Tunnels Obtained by the Substitution of an 
Equivalent Homogeneous Boundary for the Discrete Slots. NACA RM L53E07b, 

1953. 

10. Coleman, Robert P.; Feingold, Arnold M.; and Stempin, Carl W.: Evaluation of the 

Induced-Velocity Field of an Idealized Helicopter Rotor. NACA WR L-126, 1945. 
(Formerly NACA ARR L5E10.) 

11. Castles, Walter, Jr.; and De Leeuw, Jacob Henri: The Normal Component of the 

Induced Velocity in the Vicinity of a Lifting Rotor and Some Examples of Its 
Application. NACA Rep. 1184, 1954. (Supersedes NACA TN 2912.) 


84 



12. Heyson, Harry H.: Wind-Tunnel Wall Effects at Extreme Force Coefficients. Paper 

presented at International Congress of Subsonic Aeronautics (New York), Apr. 1967. 

13. Matthews, Clarence W.: Theoretical Study of the Tunnel -Boundary Lift Interference 

Due to Slotted Walls in the Presence of the Trailing-Vortex System of a Lifting 
Model. NACA Rep. 1221, 1955. (Supersedes NACA RM L53A26.) 

* 14. Glauert, H.: Wind-Tunnel Interference on Wings, Bodies and Airscrews. 

R. & M. No. 1566, British A.R.C., 1933. 

5 15. Staff of the Bateman Manuscript Project, compiler: Tables of Integral Transforms. 
Vol. I. McGraw-Hill Book Co., Inc., 1954. 

16. Lamb, Horace: Hydrodynamics. Sixth ed., Dover Pub., 1945. 

17. Whittaker, E. T.; and Watson, G. N.: A Course of Modern Analysis. Cambridge 

Univ. Press, 1962. 


NASA-Langley, 1968 11 


L-5233 


85 


